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ST-YEAR FELLOWS 


Extraction of topographic networks from digital 

elevation data 


Lawrence E. Band 

Hunter College 
City University of New York 
Assistant Professor of Geology and Geography 


Many projects in hydrology, ecology and related discip- 
lines require the extraction and organization of environmen- 
tal data bases by watersheds. The Extraterrestrial Research 
Division at NASA- Ames is currently engaged in research 
assessing forest productivity and associated nitrogen flux 
i.n a number of conifer ecosystems. The work is partially in 
cooperation with the National Park Service as a base line 
study of acid precipitation in the Sierra Nevada. As such, 
the study involves the extraction and integration of a 
number of data planes describing the terrain, soils, lithol- 
ogy, vegetation cover and structure, and microclimate of the 
region. 

My work this summer has involved the development of 
automated techniques to extract topographic networks (stream 
canyons and ridge lines) for use as a landscape skeleton to 
organize and integrate these data sets into an efficient 
geographical information system. The techniques are appli- 
cable to a wide range of environmental problems outside of 
the present application. The software is written in both 
Fortran and C, and should be portable to a number of dif- 
ferent computer environments with minimal modification. The 
data input consists of a standard raster elevation file such 
as the digital elevation model (DEM) produced by the USGS or 
the digital terrain model (DTM) produced by the DMA. Output 
consists of the set of connected drainage divides and stream 
lines such that the drainage divides completely partition 
the image plane into drainage polygons, which are uniquely 
associated to a link in the stream network structure. 

The topographic networks are first unambiguously 
defined in a geomorphic framework and then processes in a 
multi-step procedure which first nominates potential stream 
and ridge lines, and then successively grows and refines the 
features into geomorphically reasonable, connected graph 
structures representing fully complementary drainage channel 
and drainage divide networks (figure 1). The stream network 
is a directed tree graph, rooted at an outlet node. Its 
node set intersects the node set of the drainage divide 
graph at all stream junctions. The specific topologic pro- 
perties and relations between the two graph structures 
allows the construction of a powerful set of topological and 
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geometrical manipulations of the graphs and of the data sets 
they organize. 

a <\<l 

Figure 2 shows sample output of the drainage ^divide 
networks extracted from a DEM of mountainous terrain in Cal- 
ifornia. Field checking with Ames personnel revealed the 
techniques to be particularly effective in rugged terrain, 
with few significant errors occuring. Future work to be 
pursued with colleagues at Ames include extension of the 
techniques to regions with less rugged topography, improve- 
ment of the relational/spatial data base management routines 
and application to specific current and proposed NASA pro- 
grams . 
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Requirement for Torque Transducer 
in Flight Simulators 

Rosa Maria Dishman 
U.S. International University 
Associate Professor of 
Eng ineering/Sci ence 


Each NASA-ADFRF 3 axis flight simulator requires 3 DC torque 
motors* The motors provide force for the control stick (pitch 
axis and roll axis) and the rudder pedal (yaw axis)* 

New motors have been purchased which have many advantages 
over the original model# including higher torque* lower power 
requirements* and excellent resistance to demagnetization* 
However* the new motors exhibit higher cogging and magnetic 
friction* These effects are very noticeable to the pilot in the 
low force region. 

My summer project was to find a suitable torque transducer to 
be permanently mounted on each motor to correct the problem* The 
torque transducer output will be used to electronically 
compensate for cogging/friction. 

Physical configuration constraints were taken into account 
to minimize the redesign of the existing mechanical system* The 
position of the motor will remain fixed and the position of the 
pulley with respect to the motor will be the same. The pulley 
shaft will very likely be redesigned* The small space limits the 
size of the torque transducer to 4.0" long and 3*0" in diameter. 

Some of the torque transducer specifications which must be 
met are: 

I* The Torque Range 

A* Maximum torque with correct output is 400 in. -lb. 

B* The minimum readable torque output is .5 in* -lb* 

C. A maximum overload torque exists of 1200 in* -lb. 

II* Direction of Torque 

A. The torque is applied in both directions and reverses 
direction frequently. 

B. The maximum peak-to-peak rotational angle is 315 degrees. 
Ill* Accuracy 

A. Linearity: 3% of reading* 

B. Repeatability: 1% of reading. 

C. Stability of the zero point is essential (.1% of full 
scale) . 

D* Good accuracy is needed in the region near zero. 


IV, Other Specification# 

A, Recalibration should be infrequent, 

B, The transducer response should remain essentially flat to 
at least SO Herts, 

The torque transducer is expected to use the strain gage 
bridge t echnique with a gage factor of 2 to 4 mV/V, The input 
voltage will be at least 10 volts, 

After reseaching the problem, it has been determined that 
torque transducer technology does exist to meet the desired 
specifications, Several sources have been identified with 
existing commercial products that need minimal redesign, The 
estimated cost for 20 transducers including one time engineering 
development and production is *43,000, 

The next step is to get two prototypes and then to develop 
the electronics necessary for the transducer output to compensate 
for the cogging/f r i ct i on , By the end of summer 1985, torque 
transducers should be successfully electronically compensating 
for the cogging/f r i ct i on in at least 6 DC torque motors, 

Frequent discussions with Charles A, Wagner, my NASA-ADFRF 
:olleague, provided excellent guidance and technical assistance, 
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CONTROL ENHANCEMENT USING JET VANES 


Donald W. Kelton, Hlghllne College 


Purpose ; 

This study examines the effectiveness of utilizing a vane, 
deflected Into the exhaust stream of a jet engine, In Inducing 
forces and generating moments which can move or control the 
movement of aircraft. 


Areas of, and Reasons for Interest 

1. Spin recovery - There have been Instances where operational 
aircraft (fighters) have gone into flat spins and have been 
unable to recover using their existing control modes. A jet vane 
would be completely functional as long as the engine were not 
stalled out. 

2. Supermaneuverability - Operations analyses of fighter combat 
using short-range all-aspect air-to-air missiles now in 
development suggest need to quickly point and shoot, and have 
defined unusual maneuvers which will consistently defeat the 
currently favored tactic of pulling the tightest possible turn. 
These tactics involve pitching up to extreme angles of attack, 
well beyond the stall boundary, then yawing severely to dive down 
on the adversary. (Ref. 6,7,8) 

Navy studies (Ref. 2) indicate that the maneuvering enhancement 
from utilizing jet vanes can provide a 4 to 1 margin in M time-on" 
(i.e., in position to sucessfully fire upon adversary) in 
representative fighter encounters. 

3. Tailless aircraft - Future conceptual aircraft configured 
without conventional aerodynamic tail surfaces for stability and 
control have been proposed. Possible advantages offered are 
reduced drag and control simplification. A large scale test has 
been conducted in the Langley 30' x 60' tunnel with a tailless 
aircraft and both stability and control were demonstrated up to 
angles of attack of 90 . 


Scope of Study : 

The effort this summer has been to investigate the feasibility of 
employing a jet vane for yaw control. The yaw axis was chosen 
for two reasons: (1) applicability to spin recovery and (2) 

relative ease of mounting on several types of current fighters. 
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The study involved the following main avenues: 

1. Obtain tailpipe exit jet stream data from a representative 
jet engine (used 6.E. F-404 turbofan, engine of the F-18) to 
determine unit forces which can be produced on the vane. Examine 
various Mach No. /altitude flight conditions. 

2. Examine thermal effects - Determine temperature of gas on hot 
side. Select materials capable of strength retention at elevated 
temperatures. Obtain computer analysis of time- temperature- 
thickness variation for various materials, various flight 
conditions . 


3. Determine vane size through extracting data on airplane 
characteristics and derivatives and engine parameters from flight 
simulator runs, wherein a flat spin was established and recovery 
initiated . 


4. Examine configuration parameters for design and installation 
of jet vane. 

a. Possible Installations when jets are widely spaced or 
closely placed. 

b. Requirements for two axis (pitch as well as yaw) system. 

c. Possible actuation methods. 
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Efficient Analytic Commutation of Dor 1 vat Ives 
of tha Potential Energy Surfaea for Small 
Molecules or Holacular Clusters 


Harry King 


State Uni vara Ity of Haw York at Buffalo 
Professor of Ch*«tstry 


Tha anargy of a molecule or of a aolacular cluatar Is daacrlbad by Its 
potantlal anargy aurfaca (PES), a quantity which la of fundamental Importance 
for a wlda varlaty of chemical problama auch as for tha pradlctlon of chamlcal 
reactions and optical propartlas of aolaculas In shock fronts craatad by a 
raantry vehicle. Typical molecule vibrations corraspond to small amplltuda 
vibrations In tha vicinity of 1 minimum on the PES* and infrarad spectra 
corraspond to transitions batwaan quantised states of such motion* For this 
application ona Is particularly Interested In tha curvature of tha surfaea at 
the bottom of a potential well, l.e. second derivatives of tha anargy with 
respect to displacements of nuclei from their equilibrium positions* A chemical 
reaction corresponds to motion from one local minimum on a PES tc another. That 
rate of the chemical reaction Is profoundly affected by tha location and height 
of tha transition stata, l.e. the saddle point connecting these two local 
minima. Tha computational problem of finding such transition states Is 
significantly reduced If ona can computa not only tha energy for a given 
geometrical arrangement of nuclei but also derivatives of that energy. 

For a molecule with N nuclei the PES Is a function of 3N Independent 
variables , 

E - E(xl ,yi ,zl ,x2,y2 f . , . ,yR,*N) 


where (xK f yK,*K) are the three cartesian coordinates of the Kth nucleus. 
At any point on the surface, l.e. for a specified set of 3N coordinates, the 
value of the energy ,E, Is given by solving the electronic Schroedlnger 
equation. This Itself Is a formidable computational problem that has been 
addressed In thousands of articles In the scientific literature. The equations 
are never solved exactly, but from theory and experience one knows pretty v/el 1 
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how much accuracy to axp#ct fro# varloua approximate solutions* Th# 
computational chemistry group at Ames if on# of th# world* a loading c#nt#ra for 
thi* work* and th#y hav# d#v«lop#d stato-of-th#~art computer cod#a on th# Cray 
and Cyb#r supercomputers for #olacular #l#ctronic structur# calculations. With 
th# exception of on# such cod# callad GRADSCF chese programs all compute th# PES 
point by point. At th# baginning of my project th# GRADSCF cod#a w#r# capabl# 
of computing first ’erlvatives of th# #n#rgy within th# self-con* 1 stent-f l#ld 
approximation to th# Schro#ding«r aquation* i.a. 

G ■ dE/dx 


My project is to axtand this to sacond darivativas. Th# technology for 
this has* in fact, already baan davalopad at other institutions* My own 
research group at Buffalo has devoted several man years of effort to th# 
development of clos#1y related codes, but for various reasons we undertook th# 
task of writing an entirely new computer program* It Is structured to b# #or# 
compatible with existing codes at Ames, better suited to th# vector processing 
capabilities of the Cray computer, and organised with a view toward use with 
othar more accurate approximations to the true solution of th# Schroedinger 
equation such as the MCSCF and Cl formalisms. 

The quantum mechanical expression for th# electronic energy has b#«n 
d i f f #rent i ated analytically and cast in a conv#n1ant mathamatical form for 
computation. This reduces to an arduous exercise in calculus* A few thousand 
lines of new code were written during the summer for evaluating the resulting 
formulas* Just as this first summer draws to a close the computer program has 
been completed and tested and found to give results In agreement with previous 
results obtained at Ames from a series of polnt-by-polnt computations. Thus we 
are confident that our mathematics is correct. Much remains to be done to 
improve computational afficiancy and to enable the program to handle larger 
molecules. The program will then be linked up with other programs to comput# 
Infrared spectral properties and to search for transition states. 

Although I hav# d#vot#d #uch of #y professional lifetime to computat ional 
problams In chemistry, I had no prav tous #xp#rienc# th# Cray or with vector 
machines in general. So this has baen an entirely new aspect for me. I have 
worked ciosely with Dr Andrew Komornlckl throughout th# summer and h# has 
accepted th# greatest part of th# burden of educating m# in all this, but Drs. 
David Cooper, Charles Bauschl 1ck#r , Richard Jaff#, and others hav# also been 
very helpfull not only with computational mattars but also In giving m# soma 
insight into the general methods and applications being pursued within th# 
computational chemistry group at Amas. 

I hop# to continue this work n#xt year. Just working out th# calculus and 
producing a working computer program has required an intensive effort which left’ 
little tim# for thinking about chamlstry. Like a boy with a new knif# looking 
about for a good pieca of carving wood, I look forward to seeing th# program 
applied to th# solution of an actual chemical problam. 


Fault Tolerant Features of RAMPS (Redundant Asynchronous 
Microprocessor System) 


Stella Lawrence 
Professor 

Bronx Community College of the 
City University of New York 


RAMPS is the architectural structure of the hardware for UFTCS 
(Ultrareliable Fault Tolerant Control System) . In fault tolerant 
systems, decisions depend upon a redundant logic majority vote. The 
basic criteria are that the redundant system elements must be 
independent and the redundantly generated outputs must be unambiguous. 

RAMPS consists of n, autonomous, asynchronous computers in 
parallel. The asynchronous system elements have a knowledge of real 
«•*■»» and therefore generate outputs that can be correlated on the basis 
of amplitude and time. 

Systems with free integrators have control algorithms that are not 
asymptotically stable. The independent, parallel computers are subject 
to digitizing delays, and for a unit step input the two computers showed 
a constant discrepancy in value. Simulation demonstrated that if a lag 
element (filter) were substituted for the integrator the two outputs 
converged in time. Furthermore, since the Laplace transfom for the lag 
is and the Laplace tranform for the integrator is 1/s, the lag 

approaches the integrator as w approaches zero. 

My project was to solve the problem analytically. X bad to 
determine the correct representation of the control loop, which involved 
an analog input, A/D conversion, sampling, sample and hold circuits, and 
D/A conversion, obtain the system outputs y(t) and the computer outputs 
f(k), determine whether they converged and if so, how rapidly. 

RAMPS essentially involves a hybrid, sampled-data negative feedback 
system. The system is a hybrid system because the analog input is 
sampled and undergoes A/D conversion. The computer represents a 
discrete system and the plant (the helicopter) a continuous system. The 
system was modelled mathematically via the z-transform. (The 
s-transfoxm converts difference equations to algebraic equations). With 
the z-transform the response values for the continuous system are given 
only at the sampling instants, 0,T,2T, etc. the second computer may be 
considered to be delayed in time: jj^T, T 4^T, 2T + £T, etc. To obtain 
the response within the sampling interval t = (between 0 and T) , for 
the second computer, the advanced z-transform was used. 

For the first computer: 

Y(z)/R(Z) - G c (z)G(z)/ £l + G c (z)G(z)]J , and 
E(z) * Y(z)/G(z) . 
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For the delayed computers 

?(z)/R(z) * G r (z)G*(z)/ | h<J c (z)G(z)J , where G*(z) represents 
toe advanced z-t/ransform. 

The response to a unit step was found to be: 
y(t) ■ 1 - ae~ bt sin(ct +^>), 

* 

whore G (s) * 1/s, (Integrator) . 
c 

The response for the second computer (advanced z-transform) was 
found to be: 

jit) - e”^ T - ae "” bt e“^ T sin(ct + <^), as 

„ at 

T •■fcO, e *• »-^ i, and the two responses coincide. 

* 

The computer output was found to be 1+ me -bt sln(ct + ^) . 

Further work 1 a required; 

1) to determine the response of the system to a unit step, 

' (and to other inputs), for a filter (or lag function): 

G c (s) - l/(s+w), 

2) for further theoretical control loop analysis of the 
system. 

In addition, the inverse z-transform is usually obtained from the 
table of inverse z-transforms via partial fraction expansion. Computer 
programs should be developed to calculate the value of the resulting 
infinite series when the inverse z-transform is obtained by means of the 
long division of two polynomials. The LCAP program (Linear Control 
Analysis Program) should also be applied to help with the polynomial 
subroutines involved in determining the inverse z-transform. 

I found the work on the project stimulating and rewarding. I 
greatly enjoyed working with Dr. Larry Webster under the direction of 
Branch Chief Dr. John Bull. Dr. Jim Howard was also very helpful. 
Everyone was very kind and cooperative including Peggy Brown, Elaine 
Larsen, and the Library Staff. 
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Development of Finite Element 
Numerical Methods and Its Applications 
to the X-29A Aircraft 


By John B. Mcklnney, P.E. 

Cal Poly Pomona 

Associate Professor of Engineering Technology 


An Aerospace Structural Engineer is faced with the difficult task of 
designing a lightweight structure that can withstand the severe static and 
dynamic loading conditions imposed on such structures. In particular the 
X-29A is an advanced concept prototype fighter of the 1990 's, which 
utilizes forward swept wings and composite construction. 

The Aerostructures group of the Ames Dryden Flight Research Facility 
to which I am assigned is responsible for the structural analysis of the 
X-29A, specifically, I worked on the structural dynamics of the forward 
swept wing. 

In general, because aircraft structures are so complex and the 
accuracy of any analysis is paramount, the Finite Element Method is used to 
do the structural analysis. In the static case, we are concerned with 
internal forces, stresses, and deflections. In the dynamic case we look at 
the natural frequency and mode shapes of the structure which become an 
essential part of the flutter analysis. 

In order to obtain such results for the forward swept wing of the 
X-29A, we modeled the composite upper and lower wing surfaces and Internal 
structure as being structurally equivalent to a plate element incorporating 
beam elements around its edges. In all, it took 981 elements to model the 
semi -span of the X-29A wing. 

In order to solve for the free vibration analysis of the wing, the 
STARS program was utilized. In this program the well known differential 
equation of motion for a damped structure, KU + CU + MU * 0, is reduced in- 
matrix form to the solution of the eigenproblem, (K - A^M)y = 0, In which 
M and K are the mass and stiffness matrices, A is the eigenvalue which Is 
related to the natural frequency by u) = fjC and y being the eigenvector 
which gives us the mode shapes. 

Because the wing and for that matter, most all aerospace structures 
have rather large associated mass and stiffness matrices, the solution of 
the eigenproblem can be rather time consuming in terms of computer time. 


In order to reduce the solution time, we have come up with an Improved 
numerical technique for finding eigenvalues and associated eigenvectors. 

In the old method the roots (eigenvalues) were first Isolated by a repeated 
bisection procedure and Sturm sequence count after which the roots and 
vectors are determined by an Inverse Iteration method. In the new method 
the roots are Isolated as before and at the same time the roots and vectors 
are determined simultaneously. As a result, the solution times have been 
reduced by 30* to 50*. 

The attached figure shows both the undeformed structure (solid lines) 
and the shape of the first mode (dotted lines) for the x-29a. Also 
Included therein are the natural frequencies of the first six modes, up to 
15 modes have been successfully found using this new method with a 
substantial savings In time. 
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Characterization of the DNA from the Cyanelles of Glaucocystis nostochinearum 


0. Tacheeni Scott 
Northern Arizona University 
Assistant Professor of Biology 

My proposed project for the summer entailed the determination of the genomic 
size of the cyanelles of the enigmatic "algae* 1 Glaucocystis . For a century, Glaucocystis 
successfully defied a taxonomic assignment within established algal system?tics. 

It has been assigned at various times to the Schizophyceae ("blue-green algae") , 
red algae, dinoflagellates , green algae, and Glaucocystideae, a sub-class created 
for the taxonomic recognition of algae similar to Glaucocystis . A German researcher 
proposed in 1923 that Glaucocystis represented an endosymbiotic association involving 
intracellular cyanobacteria and a colorless green algal host. Few accepted Geitler's 
idea and the taxonomic position of Glaucocystis remained unsettled until the present 
year. I and my former Ph.D. advisors published the first unequivocal argument that 
Glaucocystis is indeed a unique association of cyanobacteria ( * cyanelles) and a 
colorless plant cell, probably of the green algal group (eukaryotic). 

If the endosymbiotic argument is to be substantiated beyond the remnant bacterial 
cell wall argument that I presented (Arch. Microbiol., late summer, 1984), I must be 
able to demonstrate additional arguments that the endosymbiont is truly a former 
cyanobacterium which is still identifiable as a prokaryote. Therefore, I came to 
NASA-Ames with the idea that I wanted to isolate the DNA from the cyanelle and to 
attempt to determine its molecular weight. My reason for wanting to do this was 
simple: Is the genome 10® or 10° daltons? If the molecular weight is closer to 

the lower value, then the endosymbiont could be argued to have lost some of its 
DNA since assuming the intracellular position. If the the value is closer to the 
higher number, then my argument that the cyanelles is STILL identifiable, based on 
DNA content, as a bonafide photosynthetic prokaryote. 

The project was slow in getting off the ground because my cultures were not in 
the best physiological condition in late June. I had not worked with Glaucocystis 
in over 4 years and I had decided to get around its slow growth by attempting to grow 
it in several novel fashions including bubbler cultures, shaker culture, agar cultures, 
and "spinner-bottle" cultures. The association is slow-growing and the summer simply 
slipped by before I had enough healthy cells to work with. The host cell must be 
fractioned with a glass tissue homogenizer in order to free the cyanelles for isolation 
work. 


I used standard molecular biology methods to isolate the cyanelle DNA. I was 
not however able to perform agarose gel eletrophoresis for the purpose of getting a 
rough first order idea of the molecular weight. I also attempted to isolate a plasmid 
fraction. The fractions are in buffer at this time and will have to be determined 
when I return to my institution. 

The summer was not lost however as I gained insight in many valuable discussion . 
with NASA personnel and researchers at other universities who presented seminars at 
both Ames and San Jose State University where the NASA-sponsored Planetary Biology 
course was taught this Summer, t was able to present a formal seminar at NASA in 
July. I received a good response and was encouraged by NASA's interest in Glauco- 
cystis . 


Vorticity Interaction 
The Measurement of Vortex Strength 

H. Jerome Shafer 

Technion-Israel Institute of Techlology 


The blade-vortex interaction problem, being of a most complex nature, can 
probably be best understood by gaining an insight into the basic phenomena of 
which the problem is composed. The Stanford Aero/Astro Dept, study of vorticity 
interaction* proposes to study single vortices produced by a shock tube in an effort 
to gain a fundamental understanding of vortices per se, especially in the compressi- 
ble flow regime and then to study how these vortices interact with solid boundaries. 

The shock wave leaving an open ended shock tube is followed by a flow of the 
gas traversed by the shock wave in traveling down the tube. In the current inves- 
tigation both the shock and the flow are confined in one direction normal to the 
flow direction by transparent plastic walls which as act as extensions of one set 
of walls of the shock tube. The shock and the gas are free to expand and flow 
in the other direction normal to the axis of the shock tube. A pair of symmetrical 
vortices are generated at the abrupt opening of the shock tube by the outflowing 
gas. A picture of the density at each point in the field is obtained by holographic 
interferometry. In this technique a double exposure hologram is made, first a 
tare exposure and then an exposure of the phenomenon. A photograph of the re- 
constructed hologram gives a picture of constant density contours over the whole 
field, identical with the isopycnal technique of conventional interferometry. 

One of the basic quantities needed to be determined in this study is the 
strength of the generated vortex. This can be accomplished by measuring the 
tangential velocity of the fluid in the vortex and its radial distribution. Because 
of the transient nature of the phenomena in the current investigation only velocity 
measuring methods of high frequency response can be used. 

Two possible methods of measuring the transential velocity are the Laser 
Doppler Anemometer (LDA) and the Hot Wire Anemometer (HWA). Both of these 
methods which measure at a "point" will require a previous knowledge of the vortex 
trajectory so that the HWA probe or the LDA measuring volume can be placed prior 
to shock tube firing. There is no possibility of optimizing adjustment during 
the actual experiment. The times involved in individual experiments are of the 
order of tenths of a millisecond. 

The LDA has certain inherent advantages: it is nonintrusive and as such does 

not interfere with phenomena under study; it has no elements that can be damaged 
by the flow; it measures velocity directly. Seeding seems to present no problems. 

It has been shown*t'hat dust particles of 0 . lqm to 0 . lynv require a time of about 
4 usee to be accelerated to full flow velocity after the passage of a shock wave 
of M =1.4 (Vo^200 m/sec). However, two problems tend to mitigate against the 
use of LDA: T) its use simultaneously with the holographic interferometry would 

be complicated, (2) it would require equipment not readily available in the 
shock tube laboratory. 

H.W.A., while being readily available in the shock tube laboratory has certain 
inherent disadvantages: 1) it requires an intruding sensor which may interfere 


*Bershader, D., Progress Report & Proposal for Continuation of a Vorticity Interaction 
Study; Dept, of Aero/Astro, Stanford University, July 1984. 

irk 

von Stein, H.O. and Pfiefer, H.J., "Investigation of the Velocity Relaxation 
of Micro-Sized Particles in Shock Waves Using Laser Radiation"; Applied Optics, 

Vol . 11 #2, Feb. 1972, pp. 305-307. 


H. Jerome Shafer 


with the flow, 2) the delicate wire on the probe can be easily damaged, 3) it 
measures mass flow rather than velocity directly, 4) it may present problems In 
determining relative flow direction especially In the transonic regime. The 
overcoming of these diffidulties is a major part of the effort to measure the 
radial distribution of the tangential velocity of the fluid in the vortices. 

In the shock tube experiments, as mentioned above, two symmetrical vortices 
are formed. One will be used as a reference and the other will be "experimented" 
with. A splitter plate placed mid-way between the vortices along the axis of the 
shock tube will prevent interaction between the vortices even when one is being 
"experimented" with. In a preliminary experiment a H.W. probe were placed at 
some predetermined position relative to the trajectory at one of the vortices and 
an interferogram of the complete flow field was made at such a time that interaction 
between probe and vortex has taken place. The two portions of the interferogram - 
the undisturbed vortex and "experimented" vortex would be compared to determine 
the extent of the interaction. The H.W. A. probe has been examined to determine how 
well it holds up mechanically under the forces imposed by the shock wave and flow. 
There seems to be no problem in that respect. 

It is anticipated that the main mechanical problem will be wire breakage. 

This can probably be overcome by the use of epoxy-backed sensors such as have been 
developed at Ames. Experience there has shown that the epoxy backing may 
not affect the operation or characteristics of the wire. *Robinson, et al reports 
"the epoxy-backed wire, which is useful in high speed flows because of its improved 
strength gave very similar results to the free "hot-wire". The added rigidity also 
enables a more precise determination of relative angle between flow and wire. 
Arrangements have been made to have such epoxy-backed hot-wires made at Ames. 

One of the main problems in using HWA to measure velocity in compressible flow 
is the fact that the HWA responds to mass flow and temperature. In the incompressible 
flow regime changes in flow velocity are not accompanied by changes in density or 
temperature and hence the constant value of these quantities can be combined into a 
calibration constant. The determination of velocity then follows simply from a 
direct calibration of bridge current or voltage at a function of velocity. In the 
compressible flow regime changes in velocity are accompanied by changes in both 
density and temperature. Special calibration and operating procedures such as 
repeated measurements at several overheated ratios must be followed to allow the 
separation of the several variables. 

Horstman and Rose have shown that sensitivity of the hot-wire to ambient 
temperature variations is reduced by operating the wires at as high an over- 
heat ratio as possible. An additional advantage to operating at high over-heat 
ratios is that the density sensitivity and velocity sensitivity ratio is con- 
stant. In the proposed study this is an advantage since density can be determined 


*Robinson, S.K., Seegmiller, H.L. & Kussog, M.E.; Hot-Wire & Laser Doppler 
Anemometer Measurements in a Supersonic Boundary Layer, AIAA Paper AIAA 83-1 723-, 
Danvers, Massachusetts, July 1983. 

** Horstman, C.C. & Rose, W.D., Hot-Wire Anemometry in Transonic Flow, NASA 
TM.X-62495, December 1975. 
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independently from the i nterferograms which allow a direct determination of 
the density at every point in the field at the Instant at which the second 
exposure of the hologram is made. A series of Independent i nterferograms 
made at time lags of several tens of microseconds apart will enable a 
time history of the complete density field to be made. Densities so determined 
can be used to obtain velocity from the HWA mass flow measurements. 

The technique proposed to measure the tangential velocities is that of 
the X probe. The probe is set parallel to the shock tube so that the X is symmetrically 
oriented with reference to the horizontal velocity vector. By subtracting the 
signal of one wire from the other the normal velocity component of the complete 
flow will be determined, i.e., tangential velocity of the vortex plus or minus 
the normal component of the mean flow velocity. The tangential vortex velocity 
in the direction normal to the shock tube will change direction after transver- 
sing the vortex core while the normal velocity of the general flow will continue 
in the same direction. Assuming a symmetrical vortex the vortex tangential 
velocity can then be determined by subtracting measured combined velocity 
symmetrically oriented with respect to the vortex core. 

The effect of the angle <t> between the flow velocity vector and the normal 
to a hot wi re*is given by cosnty, where m is a function of Mach number M , 

Sinits to Muck show that m = 0.5 for M< 0.5; m = 1.2 for M>2. In the range 
0. 8<M<1 . 5 m is a rapidly varying function of M . Since the whole technique 
of decomposing velocity vectors with the HWA depends on cos m <f> { the Mach number 
dependence of rn will have to be investigated further and appropriate calibra- 
tion procedures developed. 

Despite the many problems in using the HWA to and in determining vortex 
strength it is felt that measurements can be made so as to further evaluate 
the magnitude of the various parameters with at least some moderate precision. 




Smits, A.J. and Muck, K.L., Constant Temperature Hot Wire Practice in Supersonic 
Flows, Part 2: The Inclined Wire; Experiments in Fluids 2, 1984, pp. 1-9 


A Quadratlcally Convergent Scheme with 
a Large Degree of Parallelism for 
The Solution of Linear Matrix Equations 


Dan Tylavsky 

Arizona State University 
Assistant Professor of Engineering 


Solution of linear matrix equations by parallel methods Is necessary If 
large problems which are of interest to NASA are to be solved. Parallel 
algorithms to solve these problems are, therefore, necessary. Methods 
which show superior performance are those which minimize precedence 
relationships. A good example of such a method is the quadratlcally 
convergent algorithm originally proposed by Hotelling [1943]. This method 
* Is, however, limited by the restriction that the eigenvalues of the 
Iteration matrix have moduli which are less than 1.0. A quadratlcally 
convergent algorithm which has been investigated this summer Is known as 
the Approximate Multiple Factorization (AMF ) method. This method is 
superior to Hotelling's method in three ways. 

1. The method converges to the Identity matrix rather than the dense 
Inverse matrix and hence, preserves sparsity in the Iteration 
procedure. 

/ 

2. The region of convergence of the method is substantially larger than 

Hotelling's method, 

♦ 

3. The method requires one matrix multiplication per Iteration rather 
than two. 

The AMF method is competitive with Hotelling's method in two ways. 

1. The method is quadratically convergent. 

2. Few precedence relationships exists since the method represents the 
inverse as a product of sparse matrix factors. 


The first step in the AMF method is to order the coefficient matrix so that 
the largest matrix elements are located on or near the diagonal. The next 
step is to approximate the coefficient matrix as a block diagonal matrix by 
ignoring the off-diagonal elements. The inverse of this approximate matrix 
is then used as a preconditioner and will generate the sparsity pattern 
shown in (1). 
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Notice that some fill-in will occur at this point In the algorithm. If the 
system of equations is loosely coupled then an approximate factorization of 
this matrix is, 
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where the order of factorization is chosen arbitrarily. The effect of 
choosing this order in various ways is a subject of further research. 
Taking inverse of both sides of the (2) while ignoring the order of the 
factorization gives 
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Premulti plying the' preconditioned coefficient matrix by the approximate 
inverse results in the following matrix. 


A(l,l) A(l,2) A(l,3) 

M M W 

A{2,1) A(2,2) A(2,3) 

M M M 

A(3,l) A(3,2) A(3,3) 


This matrix will, from past numerical experiments, be a closer 

approximation to the Identity matrix. If the procedure is repeated then 

the results at the end of the second iteration will be even closer to the 

Identity matrix. Many numerical examples have been examined during the 

summer. One example which gives results which are typical of the 

algorithms performance, uses the matrix shown in (4). The tri diagonal 

matrix shown in (4) was chosen for this example since this type of matrix 

topography Is typical of the matrices encountered In CFD type problems. 

chosen to represent a particularly difficult 

extreme, Is Indicative of the CFD trl diagonal 

diagonally dominant nor positive definite. 
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The AMF method as described above was applied to this matrix by using 2 by 
2 and 4 by 4 diagonal block preconditioners as shown In (1) and using 
approximate factors as shown In (2). The solution vector, x, was assumed 
to have all elements values equal to 1.0. This was then used to derive the 
b vector. A graph of the Euclidean error norm of the solution estimate 
versus iteration number is shown In Figure 1. This figure shows results 
for both the use of 2 by 2 and 4 by 4 diagonal block preconditioners. The 
Important results to be taken from this figure am 1) the algorithm, shows 
quadratic convergence characteristics and 2) the use 2 by 2 diagonal 
block preconditioners has a slower Initial convergence rate, but produces 
the converged solution in the same number of iterations. The second 
observation Is Important because when 2 by 2 preconditioners are used, the 
number of calculations which need to be performed are significantly less 
than when 4 by 4 preconditioners are used. 

During this summer I have had many useful discussions with Ken Stevens 
about the role of parallel processing in the problems of interest to NASA. 
Based on these discussions I have plans to continue research on this method 
when I return to ASU, and to use the method presented as a preconditioning 
technique for use with conjugate gradient methods. 

The second project in which I have participated, is the creation and 
submission of a research proposal to AFOSR for support of parallel 
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processing research on problems of Interest to NASA. Cliff Rhoades and Ken 
Stevens have been particularly helpful by Introducing me to potential 
sponsores at AFOSR and by aiding me In constructing a proposal which Is not 
only of Interest to AFOSR and me, but also of Interest to NASA. 

References 


H. Hotelling, "Some New Methods In Matrix Calculation," Ann. Math. Stat., 
14, (1943), pp. 1-34. 
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HING/VORTEX INTERACTION 


Sam Centres 
Associate Professor 
Boston University 


Unlike most NASA-ASEE Summer Faculty Fellows, I had not selected or 
been assigned a project when I arrived at Ames in June. My NASA colleague, 
Dr. Charles A. Smith of the Low Speed Aircraft Research Branch, had given 
me a list of topic areas of interest to his Branch, and I spent most of the 
summer selecting from these a subject of interest to me, and designing a 
project appropriate to my background and to the time and resources 
available for the job. This required a considerable amount of reading, and 
numerous conversations with NASA personnel, most of whom where in the Low 
Speed Aircraft Research Branch or in the Army Aeromechanics Laboratory. 

It was established at the outset that, should a suitable project be 
agreed upon, funds would be provided by NASA Tor continuing the work at 
Boston University during the forthcoming academic year. 

The project agreed upon is a computational investigation of the 
interaction of a convected line vortex with a wing of finite span. The 
line vortex is initially straight and is aligned with the leading and 
trailing edges of the wing. As 1 t passes by the wing it will be distorted, 
and it is this distortion that is the focus of the project. The situation 
is closely related to the interaction of the trailing tip vortex left 
behind by one blade of a helicopter rotor with other blades of the same 
rotor. Such helicopter blade/vortex interactions are responsible for 
diminished ride quality, and for increased vibration and noise. 

The advantage of studying the wing rather than the rotor is that it 
reduces the extent of the modifications required of the computer code which 
will be used (one developed by a colleague at Boston University) without 
altering f undamen tally any essential aspect of the interaction. The 
existing code will have to be modified to incorporate a line vortex and to 
properly model its convection. The challenge will be to devise a scheme to 
handle the distortion of the wake as i t passes by the wing and interacts 
with the wake generated by the wing. 

Computer code modifications and initial computations will be carried 
out at Boston University and if need be continued at NASA Ames. (Similar 
computers are available at both sites.) The work will be completed and 
written up during my second stay at Ames next summer. 
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Botany Facility Feasibility Study 
For Fairchild Loaaocraft Frooflyor Plat form 

Donnie host lake 

•outhwestem Colloft Chula Vista, Calif. 

Tho purpose of this raaa a rch was to examine the critical parameters control 1 inf 
the feasibility of a plant growth facility supported as a secondary payload for the 
L o a ao c raft plat fore. The present Plant trowth Units (P9U) and the proposed European 
•pace A ge n cy (EM) Botany Facility for Europe's Free Flying. Retrievable Carrier 
(EUAECA) were studied for possible technology applications and insight. A detailed 
r epor t , describing the details of the study is available upon request. 

The Loasocraft plat fore is Fairchild's Multimlssion Modular Spacecraft (MM8) 
adaptation for s pace ceevsrcialisation vent u re. Eetieated launch dates are late 1987 
and aid 19M. The seco n dary payload specifications are listed below. 

84 cubic feet 

ISM Kg 

ISM watts (eaxieue continuous) 

Closed loop systee (negligible venting) 

ft eonth service Intervals 

Cost i S.S to i.S eillion del lars/nonth 

Leaaecraft Primary Payloads were studied to evaluate iepact on heat dissipation 
capabilities of a botany facility. A typical low earth orbit was eodelsd and energy 
fluxes for varying attitudes were calculated. A "target" po w er range of 858 watts 
was proven to be feasible within a radiative surface t sapor at are of 38 to 4S°C. 

Major requirements for a closed-loop life support system for longterm plant 
growth were divided into the following categeriesi 

Lighting 

Life Support System 
Image Collection and 
Control Systems 


Total 88ft watts 

The plant growth chamber design included two separate volumes with alternating. 
18-hour day and night cycles to effect a continuous light and heat load througout the 
orbit. Estimated volume was 8.S Ft . 

The design estimate includes two banks of high temperature flourescent lamps. 
M.9 watts of thermoelectric power to handle the plant's latent heat load, and 49 
watts allocated for water separation and recovery. A finite difference computer 
simulation was employed to examine i transient dormant phase, transient heating phase, 
and power vs radiator surface temperature for steady state orbital heating. These 
data were calculated for two different radiator masses and viewing directions. 


9ft watts 
1M watts 

98 watts 
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Figure 3, included in this abstract, ie a graph of C9S watte heating dieeipation 
during traneient phaeee for both radiator* and both sun-viewing (X 1 and Z axis 
viewing direct i one* Initial tewperaturee were C9C° K f the value approached by the 

dorwant phaee radiator* following 7 or t orbits, at the end <or beginlng) of an 
orbit. 


The Z-axia radiator* approached a conetant eteady etate value of near 9.7° C. 
Thie ie the “beet case" nodal, requiring the loweet eurface tewperature to radiate 
the heating load of COS watte. Thie direction ie non sun-viewing, "teeing" only deep 
apace, larth'e infrared radiation, and Earth's albedo. 

The cun-viewing radiator* exhibit a fact tewperature riee during the sunlight 
portion of their orbite, approaching a cyclic eteady state behavior cloee to the end 
of the fourth orbit. The significant phenomeno n here is thie short orbital tine 
necessary for the surface to increase significantly. One possible conclusion is to 
design for sun-viewingi although this is a conservative approach. If the heat 
exchange is designed for the lower tewperature off-sun viewing, a very short period 
of direct sun exposure could easily overload the systew’s capacity. 

The prograw ewployed for this study entitled, "Biwplified Shuttle Payload 
Therwal Analyser 1 * <SSPTA>, was developed by Arthur D. Little, Inc., for NASA Goddard 
Space Flight Center, Novewber 1977 and revised Septewb s r 1979. The author received 
extrewely valuable direction and supervision frow Scott Has (RFD) of ARC, who wade 
possible the prograw transfer to the Life Sciences VAX systew in building 24S-A. 

The Transient Therwal Analyser prograw solves a set of non-linear energy balance 
equations by a wodifisd Sauss-Seidel technique using N ewto n Raphson iterative 
procedure. Nany different solution techniques and convergence criterion can be 
explored with this finite difference schewe. 
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Data Conprcaioa Stadies for the Numerical Aerodynamic* 

Simulator 

Nicholas Weyland 

Department of Computer Science 
Montana State University 
Bozeman, Montana 59717 - 


The Numerical Aerodynamics Simulator (NAS) Processing System Network (NPSN) is a 
proposed computation resource being designed at the Ames Research Center. It will be capa- 
ble of handling large computational fluid dynamics (CFD) problems. Files on the NPSN will be 
quite large, up to 400 MB in the initial operating configuration. With such a large amount of 
potential traffic it is clear that any reduction in the size of files could have the potential of 
increasing throughput. This summer I have researched data compression as it affects the 
NPSN and have proposed the design and testing of a file compression algorithm for compress- 
ing the result files or solution set files generated by solving CFD problems. These solution files 
are among the largest files on the NPSN. I have determined by numerical experiments that 
they are capable, using the proposed algorithm, of being compressed by 50% with a reasonable 
amount of processing overhead. This algorithm will also compress the grid files of the NPSN. 
This will have a significant effect on the throughput and hence productivity of the NPSN. 

I will not go into the details of the proposed algorithm but refer the reader to the final 
report of the author submitted August 24, 1984 to the NAS projects office. 

I also have plans to continue this work during the coming academic year at Montana 
Stste University. I have obtained tapes of a solution file of a fluid flow problem involving the 
space shuttle and plan to use these in testing the proposed algorithm. The algorithm will be 
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designed, coded, tested, and benchmarked against these flies. 

There is also the potential far farther research into the area of data c o mpr e s si on. The 
NASA-Stanford-ASEE fellowship has provided me with the opportunity to consider data 
compression problems. The field is currently very active. 


Inhibition of Motion Sickness 
Jacob Zabara 


Temple University 

Associate Professor of Physiology/Biophysics 


A question arose as to the role of inhibitory nerves in the development of 
motion sickness, and whether inhibitory nerves can prevent or suppress 
motion-induced emesis • 

The primary approach towards the analysis of the neural regulation of 
emesis, retching and nausea is almost completely based upon excitatory 
systems involving the emetic center and chemoreceptor trigger zone. Often 
this excitation is conceived as excessive stimulation of certain receptors , 
such as the vestibular system which results in motion sickness. There 
exist extensive studies 'f excitatory systems in motion sickness or emesis 
involving latency, susceptibility, ablation of brain systems to determine 
essential areas for vomiting, stimulation of vomiting center, chemoreceptor 
trigger zone and associated regions, vestibular stimulation, visual 
stimulation and the presence of excitatory transmitters in the cerebro- 
spinal fluid. This research project will present data in relation to an 
inhibitory system which apparently acts in conjunction with the excitatory 
system. Specific knowledge as to the nature of these inhibitory nerves and 
the identification of the stimulus parameters for activation of their 
receptors, could lead to successful techniques to prevent emesis, retching 
or nausea such as occurs in motion sickness. Motion sickness studies often 
involve emesis or retching as an "end point" but the excitatory system model 
of emetic regulation has been difficult to accommodate to mechanisms of 
motion sickness (Wood and Graybiel; Grt.ybiel et al.). It may be that part 
of this difficulty has its origin in the presence of an inhibitory system. 

It has been demonstrated that autogenic-feedback training can be used to 
teach human subjects to suppress their own motion sickness symptoms (Cowings 
and Toscano). A possible explanation for this increased tolerance to motion 
sickness observed with autogenic-feedback training is increased inhibitory 
influences of higher brain structures on autonomic neural centers in the 
brain stem. However, the nature of operation of this inhibition remains 
virtually unknown. A series of investigations was initiated to identify 
and determine the specific properties and actions of an inhibitory system 
in relation to emesis, retching and nausea. There are two major goals 
of this research: (1) a demonstration of the existence of identifiable, 

inhibitory nerves preventing emesis or motion sickness and (2) the 
elucidation of the mechanisms of this inhibitory action. The initial 
research upon which this project is based, was the first to present 
significant evidence for, and identification of, an inhibitory pathway 
preventing emesis. This was based on a model of emesis produced by 
stimulation of gastrointestinal vagal afferents. I propose to extend the 
experimental model to motion sickness by analysis of emesis initiated by 
activation of the vestibular apparatus and the chemoreceptor trigger zone. 
The results of the initial research model were briefly, as follows: Emesis, 

preceded by retching, can be induced consistently in the dog by 


2 




appropriate electrical stimulation of the abdominal vagal afferent. 
Stimulation is applied usually to the dorsal or ventral branches of the 
vagus at the supradiaphragmatic level. Presumably this result is the 
consequence of generating an emetic signal in abdominal afferents which 
can oridinarily produce emesis by appropriate activation of their receptors . 
However, retching and emesis have not been reported with centripetal 
stimulation of the cervical vagus (Agostini et al.; DeBurgh and Evans; 

Evans and Murray). One possible explanation is the presence of an afferent 
pathway which is inhibitory to retching and emesis. Evidence supporting 
this explanation was obtained by stimulation of the cervical vagus 
simultaneously with the abdominal vagus. Stimulation of the abdominal 
vagus alone produces emesis or retching which is absent with stimultaneous 
stimulation of both the cervical and abdominal vagus nerves. 

The experimental procedure involves the implantation of nerve electrodes in 
animals, under anesthesia (pentobarbital, 35 rag/kg IP). These nerves had 
previously been identified as emetic-inhibitory neurons . Now, the problem 
is to determine whether stimulation of these neurons will prevent or 
suppress motion sickness. The electrode configuration has been designed as 
a nerve cuff in which the electrodes are embedded in a cylinder of silicone 
rubber. The cylinder has been slit lengthwise so that the nerve can be 
inserted with the electrode array lying perpendicular to the nerve trunk. 

The electrodes are then tested with an impedance meter to determine whether 
the coupling to the nerve is functional. The relationship between 
frequency of stimulation and impedance is non-linear so that the voltage of 
stimulation may have to be changed with a change in frequency to achieve 
optimum activation of the neural population. The measured impedance is a 
function of the electrodes, nerve and connective tissue. The impedance 
can be modelled as a simple resistance - capacitance equivalent circuit 
where the resistance is about 2.2 kjxand the capacitance about 0.50,tff. 

The capacitance is primary at low frequencies (about 10 Hz) and is due 
apparently to a double layer surrounding the electrodes . The resistance 
is primary at high frequencies (about 1000 Hz) and is due apparently to the 
nerve and connective tissue. 

Therefore, this summer we have developed the basis for experimental pro- 
cedures to explore the role of inhibitory nerves in motion sickness. 





Theoretical and Experimental Studier oi the Electronic 
Spectra o& Some Bio^lavim Adsorbed on day Surface* 


Yvonne Jones Aronowitz 
Instructor 

VeAnza Community College. 

Biologically significant organic molecules have been found to be readily 
adsorbed by clays. The amount of material adsorbed can be greatly increased 
by previously saturating the clay with certain inorganic cations. This pro- 
perty has been used by many to help explain the origin of life on earth. 

My research problem for the last two summers has been to try to understand 
the nature and extent of the bonding in these bioorganic-clay systems. 
Preliminary studies were undertaken last summer to obtain the room temper- 
ature and low temperature optical spectra of riboflavin and some simple 
bionucleotides adsorbed on clays. Changes in the electronic spectra indicated 
that there was significant interaction within the system depending on the 
nature of the exchanged cation. 

As a result of the analysis of these preliminary results, three basic 
topics naturally developed for this summer's investigation. 

The first area of concentration represents a completion of last summer's 
laboratory v;ork. It was found that the optical spectra of the cation clays 
change significantly with temperature. Therefore, the use of a room temper- 
ature clay spectrum as the reference background for low temperature studies 
introduced some error into our double beam results. An unsuccessful attempt 
was made to find a suitable instrument to repeat this part of the experiment. 
The low temperature spectra of unadsorbed cation clays were taken and these 
spectra will be used to correct the previous data. In order to facilitate 
the interpretation of theoretical results, the optical spectrum of lumiflavin 
adsorbed on clays was also taken. 

The second area of concentration was to more rigorously analyze the 
line spectra collected over the last two summers. Much time was spent 
trying to find a way to digitize the data. Equipment in the laboratory 
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of Dr. Jacob Miller is being used for this purpose. A BASIC program 
was developed to determine transition frequencies (wavelength maxima) 
and transition Intensities ( peak height and band widths). This was 
done by analyzing higher derivatives of the digitized data assuming 
simple Gaussian shapes. Most of the digitizing will be completed by 
the end of summer and the software program will be further developed 
to include other line shapes and spectrum stripping capabilities. 

The final area of concentration was to calculate the nature and 
extent of the bonding non-bonding interaction by developing a theoret- 
ical model which could be used to predict the observed spectra. To 
accomplish this phase, theoretical calculations were initiated using 
a modified Iterative Extended Huckel (IEHT) program developed here 
at NASA. The program is being successfully run on a NASA VAX/VMS 
computer. Because of the limitations of the VAX storage, studies have 
been limited to simple flavin type molecules and their geometric isomers. 
During the academic year, some minor modifications will be made in the 
program in order to rurjit on the larger Cray computer here at NASA. 

When this has been accomplished, calculations will be made on the 
composite flavin-clay systems. 

Discussions with Dr. James Lawless about future collaborative 
efforts have been held and have resulted in my applying for admission 
into the Ames Associate Program. 


•A. - 


2ND-YEAR FELLOWS 




LJ 


i 

»i 


l 

\ 


Broadband (Optical, Infrared, and Radio) Observations 
of Quasars and Active Galaxies 

Thomas J. Balonek 

Williams College 

Visiting Assistant Professor of Astronomy 


The physical processes by which quasars and active galaxies 
produce their enormous amounts of energy remains an Important 
unanswered astronomical question. Quasars, the most distant and 
luminous objects In the Universe, emit electromagnetic radiation over 
a wide frequency range - from radio to gairma ray energies. The 
relationship between the emission mechanisms and emitting regions 
throughout this wide range of frequencies Is still not completely 
understood. 

To Investigate the energy production, transport, and emission 
processes which are responsible for the prodigious amounts of observed 
broadband radiation, I organized a program In 1982 to observe a sample 
of 15 quasars and active galaxies at optical, Infrared, and radio 
wavelengths. Simultaneous observations are necessary In order to 
remove any effect that the time variability of the objects would have 
on the broadband frequency dependent Intensity distribution, fine to 
two week observing sessions were conducted at nine epochs during the 
past two years at eight frequencies In the optical (UBVRI) - near 
Infrared (JHK) region of the electromagnetic spectrum, and at four 
radio wavelengths (3mm-6cm). The multi-epoch broadband observations 
will permit us to follow the evolution of outbursts, and to compare 
these results with theoretical predictions. 

In addition to measurements of the total Intensity at many 
frequencies, we also observed the polarization characteristics 
(llnaarly polarized flux density and position angle) at optical and 
near infrared frequencies. By following the time and frequency 
dependent changes In the source polarization, we can probe the 
magnetic field structure of the emitting volumes. From our results, 
we can address problems such as: Is there one or multiple emitting 

region(s) responsible for the optical -near infrared radiation; Are 
these regions physically related to the far Infrared-radio emitting 
region; Do the emitting volumes change with time; and Is the magnetic 
field structure time or spatially dependent? 

During the summer 1983, I developed a computer program to 
calibrate and reduce the optical data. Using this program during the 
past year, wc completed reduction of the entire optical frequency data 
base. I presented preliminary results from this work at the January 
1984 meeting of the American Astronomical Society. I spent the first 
half of this summer completing reduction of the radio data base, while 
a graduate student a the University of New Mexico finished reduction 
of the Infrared data. 


M. Werner and T. RoelUg (my NASA-Ames colleagues for the 
NASA-ASEE program) have made measurements of a sample of quasars at 
eight infrared-microwave frequencies at several epochs. Several 
observations were made close In time to the optical-radio observations 
which my group has obtained. During the last part of the summer, I 
have begun analysis of the combined multi -frequency data base of the 
total flux density spectrum for several quasars. 

One of the most Interesting sources which we have studied Is the 
BL Lacertae type quasar 0J287 (the 287th source In the Jth section of 
the Ohio survey of radio sources!), In January 1983, this object 
underwent its strongest outburst In a decade - at radio as well as at 
optical wavelengths (Fig. 1). Our nine consecutive nightly 
observations In January 1983 Indicate that during the rise to and 
decay from maximum Intensity the spectral shape (Intensity as a 
function of frequency) at optical wavelenghts did not change (Fig. 2). 
Within the accuracy of our observations, we cannot detect any change 
In the optical spectral shape between January and May 1983, despite a 
factor of more than three change In the total Intensity (Fig. 3). 

These results indicate that the energy distribution of radiating 
particles which produce the optical emission has remained the same 
throughout this period o^ activity. 

There Is some evidence that the spectrum Is flatter at infrared 
wavelengths than at optical, which would be consistent with trends 
seen by previous investigators (Fig. 4). Combining these tesults of 
the spectral shape with the polarization variations observed In the 
January outburst (in which the polarization position angle and 
intensity were observed to change with time and frequency; Fig. 5), 
we have begun to model this outburst in terms of a two component 
source. In this model, two physically independent regions In the 
quasar are responsible for the observed intensity and polarization 
variations. By allowing these two components to vary independently, 
we can reproduce the observed fluctuations. 

Most of n\y time at NASA-Ames the past two summers has been spent 
performing reduction of the large data base which we have acquired 
from this ambitious observing program. During the next year, I plan 
to continue collaborating with Werner and RoelUg on the analysis of 
several quasars. Speciflcaly, we will be searching for statistical 
correlations between the multitude of variables (source intensity, 
spectral shape, fractional polarization, polarization position angle) 
as functions of time, frequency, and other variables. We will also be 
studying the broadband spectra and evolution of each of the sources in 
our sample. Preliminary analysis indicates that outbursts seen at 
optical wavelengths can be seen at later times at far infrared and 
submillimeter wavelengths, indicating a common emitting volume and 
mechanism. The data obtained by Werner and Roellig will be important 
in this analysis. Further coordinating observations of this sample of 
quasars may be scheduled If the results from our present study 
continue to prove promising. 
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’QUANTITATION OF HEAVY-ION DAMAGE TO* THE MAMMALIAN BRAIN? SOME 

preliminary findings 


A, B» Cox* end L. M. Kreft** 

♦Department of Rad. ''ogy end Radiation Biology, Colorado State 
University, Fort Co,. it, CO 80523 U.S.A., **Biomedical Reeeerch 
Division, NASA Ames Research Center, Moffett Field, CA 94035 U.S.A. 


ABSTRACT 

Histological preparations of brains from rabbles end alee exposed to different doses of 
various HZE pertieles or to lou-LET photons have been subjected to preliminary quantitation 
of radiation-induced morphometric changes* Ccnputer asslated measurements of several brain 
scructurea and cell types have been made using the K0NTR0N Automated Interactive Measurement 
System (IBAS, Carl Zeiss, Inc., Thornwood, N.Y. 10594 U.S.A.). New Zealand white rabbits 
irradiated at** 6 weeks of age were euthanatized 6.5*25 months after exposure to 6C) Co gamma 
pi.otons (LETqj -**0,3 M/pm), 20 Ne particles (LET^ * 35 ± 3 kcV/^a) or ^Ar partlcleB 
(LET "90+5 keV//im) , Measurements of stained aectlons of the olfactory bulbs of those 
animals indicate that the mean size (volume) of olfactory glomeruli Is reduced In a dose- 
dependent (and perhapa an LET-dependent) manner) aa aoon as 6.5 months after Irradiation. 
Differences between mean volumes of additional structure! have been noted when histological 
preparations of control mouse brains were compared with Irradiated apeclmens. Quantitation 
of lntermedlata and late changes in nervous (and other) tissues exposed to low- and high* 

LET radiations will Improve our ability to predict late effects In tissues of astronauts and 
others exposed to the radiation hazards of the space environment, 

INTRODUCTION 

Fur several years we have been studying late effects of particulate radiations In mammalian 
tissues, including the brains, of rodents and Xagomorphs /1,2,3,4/. Our published results 
primarily arp deacriptive of radiation-induced damage / 1/ , and until recently quantitation 
has been based largely either on subjective scoring of brain damage /2 ,4/ or on labor- 
intensive measurements such as cell counts (e.g. counts of necrotic neurons and neuroglial 
cells, etc.) /I, 3/. Included among our research goals are accurate determinations of the 
relative biological effectiveness (RBE) of particulate radiations typical cr representative 
of those to which astronauts have been and will be expoaed in the space environment, Ground- 
based experiment* have been performed on animal models using the BEVALAC of the Lawrence 
Berkeley Laboratory (Unlveraity of California) as the source of high linear energy transfer 
(LET) radiationa* While descriptive and subjective evaluations of ra'diation damage to ner- 
vous tissues have been fruitful ia the past (e.g. /l,2,3,4/) f objective large scale quanti- 
tation of such damage promise* to provide additional vaulable data for calculations of RBE 
values. The KONTRON Automated Image Analysis System (IBAS) provides a means by which that 
quantitation can be achieved. 

Late radiation-induced changes in mammalian tissues have been characterized in a number of 
djtferent ways (e.g, /l, 2, 3, 4, 5, 6, 7,8/) , but only recently has it become feasible to quanti- 
lv pathological damage and morpho-physiologic alterations accurately in large numbers of 
histological specimens. The development of computer-assisted automated Image analysis sys- 
tems has given Investigators the opportunity not only to measure areas and volumes of struc- 
tures in histological preparations, but also to determine degrees of gross pathological 
changes by rapid scanning of specimens with quantitation of, e.g., optical density differen- 
ces between regions. Until this technology wea developed, only tedious end relatively in- 
accurate measurements of microscopic structures could be made, end the date obtained thereby 
tended to be sparse. Now eystess ouch as the IBAS ellou image enhancement, computer con- 
trolled scanning, large numbers of accurate area and volume measurements, etc,, to be 
attained. Over the past 6-8 years, the authors Independently have accumulated histological 
specimens of control and Irradiated rabbit and mouse brains /l, 2, 3, 4/ in anticipation of 
access to a system such as the IBAS, and preliminary data obtained ualng the system indicate 
the! radiation effects can be quantified in such tissues. 
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HTTHODOtOCt 
Rabbit Irradiation* 

New Zealand whit* rabSlta ( Orvctolagua cunlculm ) War* exposed Co single acuta doaaa of ioni- 
sing radlacloo whan thay war* ^6 weak* old. The 2,3 ci dlamatar beam aatarad cha right aye 
aad merged through the left la each caa*, A portion of the forabraVn (Including both ol- 
factory bulba) wac Irradiated aa wall. The radiation* aaployed vara 1) bU Co gamma photon* 
(LET- -~0.3 kaV//jm) » 2) 365 HaV 20 Ha Iona (Bragg plateau LET^ -35 + 3 k*V//im) and 
3) 525 HaV *<>Ar Iona (Bragg plateau LETqo - 90 ± 5 k*V//ia), Particulate radiation oxpoaurea 
ware parforaad at the BEVALAC facility of the Lawrence Berkeley Laboratory (Unlvaralty of 
California). Further detail* of the rabbit Irradiation procedure* and of the atandard hlato- 
loglcal procedure* uaed to prepare brain tlaauea for evaluation nay be found In reference* 

2 1 d 1 6 § 7 and 8. 
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House bralna were kindly contributed by Dr. R.J.M. Pry and hi* assc sites from representative 
groups that had been Irradiated at doses below 10 Gray In the course of a study on high-LET 
radiation carcinogenesis /5/. All mice were perfused via the heart with a paraformaldehyde- 
glutaraldehyde fixative In aodlua cacodylate buffer. Ten-jim thick serial sections were cut 
from brains, Including olfactory bulbs, in the coronal plane and stained with hematoxylin 
and eosln* 

Automated Interactive Image Analysis 

A Zeiss Ultraphot photomicroscope was used for all tlsaue observations and measurements, A 
television camera mounted on the microscope was connected to a video monitor and computer. 
These In turn were Interfaced with e second video screen and a second computer. By means of 
a curaor and a magnetised tablet, parte of stored images were outlined, end tbe computer 
determined the ereaa wiehln those outlines baaed on calibration factors preset by the InVea- . 
tlgatora. Area values were stored on magnetic disks for subsequent statistical evaluations. 

Measurements of Rabbit Olfactory Glomeruli 

A ssglttal section through the olfactory bulb of a given rabbit brain specimen was located on 
each slide to be evaluated. Code numbers were assigned to each set of measurements to ba 
performed on each field. All measurements were made of fields magnified 300 times, and a 
single preset scale factor used for all area determinations waa entered into the computer 
before any evaluations were made. A field including several olfactory glomeruli waa found, 
and che focused image was entered into computer #1. The image was held on the screen end 
enhanced to improve contraat. Interfacing between the enhanced Image and the magnetic tablet 
allowed each olfactory glomerular area to be outlined; outlines appeared aa overlays on 
monitor #1. When all glomeruli within a given field were outlined, their areas (In ^m 2 ) 
were determined by the computer and results were stored on magnetic dlska, When adequate 
dat a had been accumulated, information fro* each animal and from each group of animals 
(based on type of radiation received, dose, and poat-lrradlatlon time) waa evaluated using 
additional IBAS statistical evaluation programs. 

Heasurementa of Structures in Mouse Brains 


Several attempts were made to measure olfactory glomeruli in the mice, but It soon became 
obvious that the delineation of representative structures was not as distinct in the mouse 
as it was in the Labbit. Therefore, since serial section* were availsble, volumes of well- 
defined rones in the bulb were determined, e.g. outer plexlform layer, granular layer, mitral 
and inner plexlform layers. Data from these measurements will be available at a later time. 
For the moment, suffice it to say chat provocative results, that require confirmation, have 
been obtained. 

RESULTS AHD DISCUSSION 

Radiation-induced Change* in Rabbit Olfactory Glomeruli 

The olfactory bulba of rodents and lagomorpha are relatively large etructuree consisting of 
regions which can be dletingulahed in etained histological preparations. The glomeruli 
comprise a well-defined leyer of "spherical" structure* In the rebbit. Sensory input from 
the olfactory mucoaa travels via processes of neuroepithelial cells which synapse in the 
glomeruli of thm olfactory bulb with processes of tufted and mitral cells /9 /, Measurements 
of individual glomerular areas on 5- to 10-yim chick sections provide an indication of glo- 
merular volumes. Since all rabbit brain sections examined were cut from similar levels /4/, 
comparisons between experimental and control groups were made after mean glomerular areas 
were obtained for each rabbit using the IBAS. 
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Erellmlnary mtaauramenu with th« IBAS ware Midi of sectioned olfactory glomeruli fro. rab- • 
blta which had racaivad ralatlvaly hl|h dotes of radiation 6*5-25 months bafora thary ware 
authaoatiaad. Haan glomerular araaa for control animals ware graatar than thoaa found for 
all irradiated rabbits Measured (aaa Table)* Host maasursmanta ware »ade on tiaauaa obtained 
17-25 Months after Irradiation! but the mean glomerular araaa ascertained for two groups of 
argon loo irradiated animals at 6*5 or 7 Months poat. exposure- also were low, implying that *" 
shfldUge of olfactory glomeruli may be a relatively early effect from which brain tissue 
cannot recover after the doses (fluxes) employed* 


TA1LK Mean Areas of Olfactory Glomeruli from Rabbits Exposed to 
Radiations of Different Linear Energy Tranafer (LET) 


Type of 

LET a 

Number of 

Doae 

Months at Risk 

Mean Glomerular Araa 

Radiation 


Rabbits 

(Cray) 

(maan) 

(^■2 k10 3) 

Control 

— 

4 

0 

10 

8.5 

60co gamma 
photons 

0.3 . 

3 

3 

10 

20 

24 

16 

6.1 * 
3.8 ’ 

2& Nt Iona 

35 ± 3 

4 

10-15 

16 

4,8 

Iona 

90 + 5 

6 

7-9 

16 

5.2 



4 

9-12 

7 

5.5 


Tha data were analyzed statistically in two different waya. All mean 

glomerular area values for irradiated animals were algnif lcantly different from mean values 
for control rabbits baaed on t-tests. In addition, Hann-Whitney U tests were performed to 
compare median valuea of glomerular areas for irradiated animals with those for controls. 

All irradiated specimens measured yielded median glomerular area values which were aignifi- 
csntly lower than control*, The dramatic reduction in sire of olfactory glomeruli from 
irradiated rebblts could stem from one or more factors. Lose of neurons from the olfactory 
bulb is possible, especially at such long intervals after irradiation. Reduced numbers of 
neuroepithelial cells end conccmaitant diminution of olfactory input might have led to the 
shrinkage of eenaory processing structures in the forebrain, but in the case of the irradia- 
ted rabbits under discussion, a relatively small proportion of the olfactory epithelium was 
exposed to ionizing radiation in each subject, Even if no neurons were lost, their ability 
to maintain pre-irradiation synaptic relationships at an optimal level might have been dimi- 
nished by irrsdlatlon, and this in turn could have led to reduced function and partial atro- 
phy. It la unclear whether the observed effect is due to the response of the glomeruli to 
direct damage or to secondary neuronal/gllal effects. A final factor could be early, inter- 
mediate end late radiation effects on the microvasculature. Quantitation of such damage .will 
be achieved moat accurately by measuring end/or counting the appropriate vessels in Che per- 
fused specimens of the irradiated mouse brains on hand. 

In conclusion, the radiations of high LET seemed to cause earlier and more dramatic shrinkage 
of olfactory glomeruli In exposed rabbit brains than did comperabk doses of ^°C b gamma photons. 
Of course, more dose groups from all experiments must be examined using tissues from animals 
euthanatized at earlier and much later post-lrradlatlon times before final statements can be 
made about Che valuea for the relative biological effectiveness of heavy iofts pertinent to 
space biology* Nonetheless, the data shown in the Table do allow us to suggest that radia- 
tion cf facta may be measured in histological preparations of brain tissues, and that by means 
of cools such as IBAS, quantitation of radiation effects In tissues previously thought too 
difficult to study is now a realistic gorJ . 
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THE USE OF ADVANCED COMPUTER TECHNIQUES IN THE 
CARDIOVASCULAR RESEARCH LABORATORY 


BY: JOHN S. DAPONTE 
COMPUTER SCIENCE DEPARTMENT 
SOUTHERN CONNECTICUT STATE UNIVERSITY 
501 CRESENT ST. 

NEW HAVEN CT. 06515 

My assignment in the CARDIOVASCULAR RESEARCH LABORATORY at NASA/AMES this 
past summer involved the use of advanced computer techniques in three applica- 
tions. The first application was the use of a sophisticated signal processing 
software package named SIC3 to perform FOUIER ANALYSIS on data obtained from 
the COSMOS 1514 PRIMATE SPACE FLIGHT. The second application was the use of 
the BMDP, statistical software package, for data reduction of experimental re- 
sults obtained from water immersion studies conducted at NASA/AMES. The third 
application was the design, coding and implementation of a program to calcu- 
late and plot various phisiological parameters resulting from ultra sound 
crystal implanted in RHESUS MONKEYS. 

The object of the signal processing effort was to calculate VASCULAR IM- 
PEDANCE as the transfer function of carotid arterial pressure and flow. After 
investigating various aspects of SIG the following general procedure was de- 
veloped: 


1* Transfer uncalibrated pressure and flow data from the PDP 11/34, in 
the Cardiovascular Lab, to the VAX 780, in the computational center, 
via the cardiovascular network. 

2. Use SIG to read the VAX data files into the SIG data base. 

3* Apply calibration slope and intercept parameters by using the CMULTI- 
PLY and CADD SIG functions for both pressure and volume time signals* 
These factors are currently available for flow and should be avail- 
able for pressure in the near future. 

4. Window time signals with a FULL HAMMING WINDOW as supplied in SIG. 

5. Transform signals into the frequency domain using the FAST FOURIER 
TRANSFORM ( FFT) function of SIG. The windowing of step 4 will pro- 
duce a frequency spectrum in which a larger portion of the total en- 
ergy of the signal is in the main lobe with side lobes that decrease 
rapidly in energy as frequency increases. 

6. Divide the frequency spectrum of the pressure by the frequency spec- 

trum of the flow yielding the frequency spectrum of the VASCULAR IM- 
PEDANCE . / 

7. The amplitude of the VASCULAR IMPEDANCE can easily be plotted using 
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the FSQUICKPLOT routine of SIG. If a phase plot is desired then the 
FSPLOT function, with changes in various defaults, should be used. 
Typing the number 6 in response to the prompt after each SIG plot 
provides for rescaling. 


A typical SIG program for producing VASCULAR IMPEDANCE would be as fol- 
lows: 


TSREAD 

1 


E270CAP. 329 

TSREAD 

2 


E270CAF. 329 

CMULTIPLY 

2 

3 

22.389567 

CADD 

3 

3 

-0.561142 

CMULTIPLY 

3 

3 

3. 19 

CMULTIPLY 

1 

4 

M 

CADD 

4 

4 

B 

FFT 

4 

5 

0 

FFT 

3 

6 

0 

DIU 

5 

6 

7 

FSQUICKPLOT 

7 




Where E270CAP.329 and E270CAF.329 are the VAX file names for the uncali- 
brated pressure and flow respectively. The values 22.389567, -0.561142 and 
3.19 are used to calibrate and convert the flow signals, while M and B are 
calibration factors (to be supplied) for pressure. 

The objective of the statistical applications was to use analysis of var- 
iance to determine if there was a significant difference among various groups 
of monkeys before water immersion, during immersion, and after immersion. 
Groups consisted of conscious, tranquilized and anesthetized monkeys. The P2V 
routine of BMDP was identified as providing the most applicable results. This 
type of experiment is called a REPEATED-MEASURES DESIGN wth one grouping fac- 
tor and one within factor. BMDP programs were written for various dependent 
variables and results have been provided. 

A'*sample BMDP control file used to direct P2V is as follows: 

/PROBLEM TITLE IS ’SODIUM EXCRETION DATA’. 

/INPUT VARIABLES ARE 4. 

FORMAT IS ’ ( FI . 0 , 3F5.2)’. 

FILE IS NAEXCR.DAT. 

/VARIABLE NAMES ARE GROUP, NA1 , NA2 , NAi. 

/DESIGN DEPENDENT ARE 2 TO 4. 

• LEVEL IS 3. 

GROUPING IS GROUP. 

NAME IS TIME. 

/END 

The primary object of the dimension processing program is to calculate a 
large number of CARDIOVASCULAR PARAMETERS such as HEART RATE, AORTIC PRESSURE, 
LEFT VENTRICULAR PRESSURE, VENTRICULAR WALL THICKNESS, AORTIC FLOW, MAXIMUM 
and MINIMUM LEFT VENTRICULAR VOLUME using spherical and elliptical geometric 
models. Another objective was to produce PV LOOPS. Thus the; dimension pro- 
cessing program outputs a file that contains intructions and data that can 
plot PV LOOPS using the TELLAGRAF software package. This program also con- 
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can be compared. expended during the different CARDIAC CYCLES 

» system flow chert of the dimension processing program £ollovs . 






SIGNAL DETECTION IN THE 

SEARCH FOR EXTRATERRESTRIAL INTELLIGENCE (SETI) 


Stanley R. Deans 
University of South Florida 
Professor of Physics 

Scientists at NASA-Ames Research Center, Stanford University, and the 
Jet Propulsion Laboratory are carrying out a five-year research and development 
program to develop powerful instrumentation to attempt detection of signals 
generated by other intelligent beings, if such signals exist. Hy contribution 
to this effort is in the software development for the signal detection system 
ultimately to be used for extraterrestrial intelligence signal searches. Hy 
work on this project can be divided into three related topics. These topics 
and associated progress reports are listed below. 

A. Development of new signal detection algorithms 

Adaptation of the Radon transform 1 to the signal detection problem was 
accomplished during the summer of 1983. Further refinement of the technique 
continued during the following academic year at the University of South 
Florida, Tampa. The remaining part of this work is closely related to the 
evaluation of detection algorithms. 

B. Development of methods for systematically evaluating detection strategies 

An important part of evaluation and subsequent optimization involves the 
rapid computation of various statistical parameters and the determination of 
their interdependence. A considerable portion of the current summer work has 
been devoted to the development of efficient computational procedures and 
related illustrations of the important interrelationships. Samples of these 
results are shown in the Figures (obtained using the DISSPLA graphics package). 

It should be emphasized that the major difficulty (and thus the major 
contribution) is the extension of existing tabulated results 2 to regions that 
are excessively tedious to compute due to overflow and/or underflow problems. 
(This is related to the very-low-probabi 1 ity searches that must be carried 
out in SETI.) It is necessary to sum hypergeometric series that contain many 
large terms, such as 3.87512 x 10 487 The trick that makes the computations 
feasible is to sum the hypergeometric series in a speciat logarithmic mode that 
produces the logarithm of the desired function rather than the function itself. 



C. Optimizing detection strategies 

Work on optimization is just beginning. Discussions with several members 
of the SET I group at Ames are underway and various approaches to constrained 
optimization are being studied. 3 * 4 Further work on this *nd other SETI related 
projects will continue during the coming year when I will be on leave from the 
University of South Florida to accept a National Research Council (NRC) Senior 
Research Associateshlp at NASA-Ames Research Center beginning October 1, 198^. 
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FIGURE CAPTIONS 

Fig. 1 Chi-square density. The probability of false alarm P. is the shaded area from 

f a 

threshold b to 00 . 

Fig. 2 Non-central chi-square density for fixed signal-to-noise ratio r and different 
degrees of freedom. 

Fig. 3 Same as Fig. 2 with different r. 

Fig. k Relationship between signal-to-noise ratio in decibels and the probability of 
missing a signal with fixed n and P. . 
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Microbiology of Hydroponlcally Grown Plants In CELSS (Control lad 
Environment Ufa Support Systems) 

Olane Dudztnskl, Ph.O. 

College of Santa Fe 

Professor of Biology 

Chairperson, Dept, of Science/Math 

Survivability of man In an extraterrestrial environment, such as a 
space station, for extended time may In part depend upon life support sys- 
tems using higher plants to supply oxygen, potable water and contamination- 
free food. Such an environment on earth can be simulated by using a CELSS 
approach In which higher plants are grown In recirculating hydroponic and/ 
or aeroponlc chambers under controlled environment conditions. Prototype 
chambers have been developed and modified over the past few years at NASA/ 

AMES In the CELSS laboratory of Robert MacE l roy/Steve Schwartzkopf . During 
this time mechanisms for controlling Indigenous microbial populations In 
lettuce and tomato monocultures have been modified and compared to deter- 
mine relative effectiveness In reducing microbial abundance and altering 
species diversity within such systems. 

Isolation and Identification of microorganisms Involved use of spread 
plate preparations of microbiological samples on Tryptfcase soy agar follow- 
ed by Inoculation of colony Isolates Into Hoffman La Rouche miniaturized and 
rapid microbiological systems. Samples were collected aseptlcally 3 times per 
week from plant troughs In recirculating nutrient systems. Nutrient samples 
were collected simultaneously with microbiological samples and were analyzed 
for micronutrients. Lettuce and tomato seeds were germinated between moisten- 
ed filter paper and planted within the growth chambers between alcohol -washed 
polyurethane plugs. Prior to planting the hydroponic nutrient delivery sys- 
tems were flushed with 20% sodium hypochlorite and rinsed with deionized 
water for a few hours. Nutrient solution was recirculated through the chamber 
until plant harvest at the end of 3~5 weeks. Control of microbial population 
density was achieved by continuous use of an on-line ultraviolet lamp emitting 
2^5 nm, or a .22 p bacteriological filter placed within the recirculating nut- 
rient system. The bacteriological filter was replaced weekly. 

. Microbial cell densities reached 10°cells per ml. to control units and 
1(T cells per ml. In UV and filter treatment units. Past problems Involving 
Iron loss due to breakdown of the Hoagland nutrient solution chelator by UV 
treatment were overcome by weekly replacement of lost Iron. Such problems 
were nonexistant in bacteriological filter treatment systems. Bacteriological 
filter units generally reached 10® cells per ml. prior to filter replacement. 
Figure 1 shows microbial cell numbers during growth of lettuce monocultures 
In a controlled hydroponic plant growth chamber. UV radiation and .22 p bacter- 
iological filter treatments reduced bacterial cell numbers by 50$. Lettuce 
plants were grown to maturity within 21 days as compared to 35 days In earlier 
experiments. Comparisons of microbial species diversity show differences be- 
tween control and UV or .22 p filter treatments, as well as between lettuce 
and tomato monoculture systems as appear In Figure 2. Future Investigations 
could be directed at studying microbial community structure within other mono- 
culture and polyculture systems and In particular at lettuce and tomato poly- 
cultures. Also of interest would be studies regarding selective seeding of 
naturally occurring and/or bio engineered microbes , as well as bacteriological 
systems modelling of indigenous and exogenous pathogenic species within such 
systems. 







BACTERIA ISOLATED FROM TOMATO MONOCULTURE BACTERIA ISOLATED FROM LETTUCE MONOCULTURE 
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CAPACITY AND CODING FOR A DIRECT-DETECTION 
OPTICAL FIBER CHANNEL 

Mark A. Herro 
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Previous research on coding for direct ’detection optical communications has primarily 
been focused on the quantum limited channel in the absence of thermal receiver noise. 
For uncoded systems, considerable analysis has been performed on models which include 
both thermal noise and an Avalanche Photo-detector (APD). In this paper we introduce 
a general non-additive Gaussian noise channel model whiJi includes the effects of both 
quantum and thermal noise. We use this model to obtain the capacity and performance of 
a direct-detection optical communication system using 2 L -ary Pulse Position Modulation 
(PPM) and an APD. The performance is evaluated as a function of the receiver threshold 
setting and the PPM alphabet size. 

We examine the performance of rate 1/2 Reed-Solomon (RS) coded PPM systems. 
The symbol size of the RS code is 2 L — 1 when 2^-ary PPM is used. The receiver operates as 
follows. The received voltage is compared to a predetermined threshold and if the voltage 
in one of the 2 L time slots in a codeword exceeds that value, a pulse is declared present; 
otherwise, no pulse is detected in that time slot. If only one time slot produces a pulse, 
the corresponding symbol is passed to the RS decoder. Two demodulating methods are 
considered for generating the RS code symbols when either no pulse or more than one pulse 


is detected in a word. In the first method, if more than one time slot in a word contains a 
pulse, one of these time slots is randomly chosen to generate the code symbol. If no pulse is 
detected in an entire word, one of the 2 L slots is picked at random to identify the symbol. 
An alternate approach is to declare an erasure except when exactly one pulse is detected 
in a word time. This information is passed on to the decoder as “side information.” The 
RS decoder then acts as an error plus erasure decoder. 

The error performance of these two schemes is compared as a function of the threshold 
setting and the number of bits, L, in a code symbol. Previous work has shown that in 
the absence of thermal noise, the best threshold (in terms of minimizing the overall error 
probability) should be about 1 photon. However, when thermal noise is present, our 
channel model predicts that the optimum threshold setting is much higher, in most cases 
greater than half the received signal level. In contrast to the noiseless channel, where less 
than one photon per information bit is required to produce decoded error probabilties on 
the order of 10 -9 , we find that even when a high quality receiver ( Z = 1500) is used, 20 
to 40 received photons per information bit are required to produce these same error rates. 
When side information from the demodulator is available, the performance is about 1 dB 
better for overall error rates ranging from 10 -7 to 10“ 11 and when 4 < L < 14. It turns 
out that the performance is exeedingly sensitive to the threshold setting. This last result 
suggests that in situations where the signal level fluctuates considerably, an alternate 
demodulation strategy (e.g. 6 max or Maximum Likelihood Decoding (MLD) ) may be 
required. An error probability comparison is made between “threshold” demodulation 
and MLD to access the penalty paid for using the simpler, but sub-optimum receiver. 

Finally, we point out that as L increases, the capacity of the channel increases but 
approaches a finite limit in contrast to the noiseless channel that ha? infinite capacity 
as L increases without bound. In addition, as the PPM alphabet size increases, more 
bandwidth is required and the overall performance becomes even more sensitive to the 
threshold setting. This suggests that a moderate value for the alphabet size may be more 
advantageous when both complexity and sensitivity to signal levels are taken into account. 
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Spt rul tna : Cultures In Shaken vs. Stable Media 


Ray A. Hill 

Morgan State University 
Associate Professor of Biology 

The on-going challenge of the CELSS (controlled environmental life support 
system) group Is to develop a reliable, functional controlled chamber In which edible 
plants can be grown aboard the spacecraft during long missions, minimizing the Initial 
payload weight, and maxlmlxlng the use of re-cyclable materials aboard the spacecraft. 
C0^,for example, can be utilized by select plants for conversion Into edible, organic 
matter, and 0 Z , on the other hand a byproduct of the process, can be utilized by 
persons aboard the spacecraft for necessary aerobic, metabolic functions. Several 
plants have been studied In detail under controlled chamber conditions to establish 
a data base on conditions necessary for maximum biomass yields, as well as control 
mechanisms for potentially harmful microorganisms that may represent endogenous con- 
taminants on the seeds and/or seedlings used to Initiate germination or planting 
(Dudzlnski, Schwartzkopf ) . Spl rul ina , a filamentous member of the cyanobactcrlal 
population, has received Increased Interest by the CELSS group at NASA as a candidate 
organism for growth’ In a controlled environment. This interest stems, In part, from 
the comprehensive report on the edible qualities of Spl rul Ina (ClferrI ,1983) , high 
protein content a chief factor. 

Becker (1978) reported the following nutritive value data of S. platensls : 

Protein Efficiency Net Protein Utilization Biological Value Digestibity 

Coefficient 


1.80 62.0 75.0 83.0 

Several studies confirm that persons living In several geographic areas have consumed 
Spl rul ina for food for many years (Ptrie, 1975; Delpeuch, F. , A. Joseph, and C. Cavelier, 
1976)". Protein content reports of Spl rul Ina grown under laboratory conditions 
have ranged from 62 to 68% of dry weight (Clement, G., C. Glddey, and R. Menzi, 19^7) . 

Current interests in laboratory culture of Spl rul Ina involves the matting that 
forms as a result of the filaments intertwining during growth. Such matting 
possibly Interferes with growth rates, I. e., minimizing protein yield, and quantitive 
research analyses. Dr. Ellen Weaver, San Jose State University has developed a defined 
growth medium for Spl rul ina , one that would minimize the matting as a result of Inter- 
twining. Dr. Richard Radmer of Martin-Marietta has also developed a defined medium 
to reduce mat formation during culture. Using the medium of Radmer at 37 ®C, Spl rul I na 
was cultured under continuous flourescent light for 5 days. One inoculum was stable, 
while the other was shaken continuously. This was done to determine if continuous 
motion would reduce or inhibit mat formation as well as formation of prevloulsly reported 
spiral and helical filaments. COj. was bubbled Into both culture flasks. Spectrophoto- 
metrlc analyses showed that stable and shaken cultures grew at the same rate, however, 
mat formation was observed In the stable cultures. Microscopic examination of the 
cultures revealed that more filamentous spirals formed In the stable cultures, whereby 
straight filaments formed more prominently in the shaken cultures. The defined medium 


seemed to reduce matting, however, the effect was more noticeable in the shaken cultures. 
The mechanism for reducing matting by using a defined medium Is thought to be related to 
salt content. Further studies are at hand to determine the relationship of salt con- 
tent to helical filament formation and spiral filament formation as compared to straight 
filament formation. The stage Is also set for futther studies on Spl rul Ina to more 
accurately determine protein content and Identify amino acid composition. 
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Separated Flows Near the Nose 
of a Body of Revolution 


Sung P. Lin 
Clarkson University 


The first project Is on the flow separation from a deforming circular 
cylinder. This problem arises in an effort to understand how three- 
dimensional flow separation tukes place near th^ nose of a body of revolution 
making an angle of attack with a uniform stream 1 . To obtain an analytical 
solution to this problem, an Iterative method of solving the unsteady 
Navier-Stokes equation was invented last summer. The method involved, in the 
last step, solutions of some integral equations which were singular. A mean 
is now found to eliminate the singularities. To check the accuracy of the 
method and computation we have determined the time history of the migration 
of primary separation points on a rigid circular cylinder. The results are 
compared with known numerical and experimental results. The comparison is 
very good as can be seen in Fig. 1. The next step is to obtain the 
separation points on a deforming cylinder. This will be followed by the 
mapping of the unsteady separated flow fields. We hope to complete this 
project by the end of next year. 

The second project is on the formation of spiral vorticies on a rotating 
body of revolution with a nose tip. The onset of vortical flows at some 
distance downstream from the nose is considered as the consequence of the 
Instability of the axisym(.'.ric basic flow about the rotating body. The basic 
flow varies rapidly in space. As a consequence, no known stability theories 
are adequate for our purpose of predicting the onset of the spiral vorticies. 
A new approach of stability analysis of non-parallel flows is formulated. An 
important realization is that the non-parallel flow instability can be solved 
completely only if it is treated as an initial value problem. The basic 
steady flow is currently being obtained by use of the Galerkin method. The 
general method is being used to generate known solutions of flows about a 
rotating sphere at small Reynolds numbers and over a disk rotating at large 
Reynolds numbers. This part of the work is expected to be completed by the 
end of this year. The basic flows about a cone and other bodies of 
revolution and their stability analysis will be carried out next year. 

Other related projects I have initiated during my stay at Ames this 
summer include the formulation of stability analysis for the flows near 
saddle points o£ separation, and the reinterpretation of three dimensional 
steady separation' 5 in the frame work of dynamical systems. 

Without the close interaction with Murray Tobak and the encouragement 
and support of Gary Chapman and Lewis Schiff, whatever progress I have made 
would not have been possible. All of the projects I studied at NASA-Ames are 
relevant to the mathematical modeling of the aerodynamic characteristics in 
flight dynamics. ’ 


n • 
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During +3G 2 Acceleration 
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Results of the analysis from last years investigation at NASA Ames 
indicated that after bedrest various physiological and demographic variables 
could not be used to predict how long human subjects could tolerate +3G Z 
acceleration. The only variable which could be used to predict +3G Z toler- 
ance after bedrest was pre bedrest +3Gz tolerance time. The question was 
then asked if pre bedrest tolerance could be predicted using the same inde- 
pendent variables. The analysis indicated that pre bedrest tolerance time 
could be predicted with very little error. The variables with tie best pre- 
dictive power included age, sex, and height. An older, Shorter male was 
found to be a beneficial profile for increased tolerance time. Variables 
with moderate predictive power included weight (heavier subjects have higher 
tolerance) , diastolic blood pressure (higher pressure was beneficial) , and 
plasma volume per kilogram of body weight (less volume was beneficial) . 
Variables that had no predictive power were those related to measures of 
fitness (yc>2 max, resting heartrate, and percent body fat) . This result 
for the fitness variables was probably the most significant finding since 
it has not yet been agreed upon as to whether increased fitness is beneficial 
or detrimental to +3G2 tolerance. Apparently other variables are much more 
important than fitness and leave little variance for fitness vari abl es to 
explain. Subjects used to create the model could be correctly classified 
into high and low tolerance groups with 90 percent accuracy. 

The value of this prediction model needs to stand up under cross-val- 
idation before its true worth is known. Using 16 new subjects who partic- 
ipated in the bedrest study for 1984 the model derived in 1983 was subjected 
to cross-validation. Results indicated almost as good a prediction under 
cross-validation as with the original data that derived the model. The or- 
iginal model predicted low and high tolerance groups for the 16 new subjects 
with 85 percent accuracy. The reduction of 5 percent is minimal and indicates 
that the original results fran last year seem to be very reliable. Additional 
subjects with varying profiles across the independent variables need to be 
tested before the final predictive power of the current model is determined. 


Although my main activity wdthin the Bicmedical Research Division was 
to determine the parameters for a model to predict +3G Z tolerance, I also 
consulted with a variety of researchers within the division. I provided 
statistical consultation for the design and analysis for current studies 
and studies to be proposed in the future. I advised programmers who were 
writing statistical programs to be used within the division and helped 
interpret results from statistical packages currently running at Ames. 


An Estimate of Aerosol Effects on the Inference of SO 2 Concentration In the 
El Chi chon Volcanic Plume 

J. L. Mergenthaler 
Assistant Professor of Meteorology 
Florida State University 
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There were no In-sltu and few remote measurements of the composition of the 
stratospheric cloud resulting from the March 28 through April 4, 1982 eruptions 
of El Chichon during the first few days of its existence. Satellite observations 
of the atmosphere's spectral reflectivity gathered with the total ozone mapping 
spectrometer (TOMS) aboard Nimbus 7 were used by Krueger(1983) to infer the SOg 
mass in the stratospheric plume on April 6 after the major eruptions of April 3 
and 4. The accuracy of this estimate, 3.3 megatons, is of great interest since 
it sets the initial conditions for the evolution of the sulfuric acid aerosol 
later observed. 

A possibly important aspect of this measurement not treated by Krueger is 
the effect of the volcanic ash and other aerosol embedded in the plume. The 
neglect of these particulates could result in erroneous S0« amounts. 

Figure 1 is a slightly adapted version of Figure 2 from Krueger's report. 

It depicts a spectral scan for five wavelength channels of the area perturbed 
by the stratospheric volcanic plume. An example is presented by Krueger in which 
this data is used to retrieve the columnar SO 2 concentration in part of the April 5 
plume, scan 92. This example is reexamined taking into account the effects of 
aerosol scattering and absorption upon the plume albedo with the purpose of 
examining the magnitude of the error in stratospheric SO? mass previously reported 
and to examine limitations on the aerosol optical properties consistant with the 
measurements shown in Figure 1. 

A three layer radiative transfer model has been developed to represent the 
reflectivity of the atmosphere containing a volcanically perturbed strata between 
44 and 22mb. In the unperturbed case only the scattering from air and ozone 
absorption need be considered above 44mb. The reflectivity of the atmosphere 
below 44 mb is assumed to be Lambertian and is computed from the off-plume part 
of the scans shown in Figure 1. In computing the reflectivity of this layer the 
climatological ozone amount has been used to adjust for the attenuation caused 
by the ozone above 44mb. The volcanically perturbed stratum is assumed to contain 
aerosol and SO 2 , as well as, ozone and air. The reflectivity of this layer and 
the underlying reflector is computed with a matrix-method doubling-adding radiative 
transfer code similar to that developed by Twomey, Jacobowitz and Howell (1966). 

It has been assumed that the reflected radiance which would have been 
measured had the plume not intervened can be approximated by interpolating the 
spectral reflectivity data across the disturbed area as indicated by the dashed 
lines in Figure 1. It is pretty clear that the reflectivity at 312.5 and 317. 5nm 
was reduced substantially by the plume; however, the reflectivity at 380. Onm 
appers to have been enhanced slightly, approximately M> figuring by the procedure 
outlined above. 

A phase function computed through the use of Mie theory with an asymmetry 
factor of .83 was used to represent the scattering pattern of rather large 
particles in the volcanically perturbed stratum. The single scattering albedo 
and optical depth of the aerosol are considered dependent variables to be 
determined from the observed reflectivities at 339. 8:*. 331 .2 and 380. Onm. 


Model results shown In Figure 2 demonstrate that the reflectivity of the 
atmosphere at 339. 8nm Is very sensitive to aerosol scattering. Figure 2 shows 
for example that a model aerosol optical depth of 1.0 of conservatively scattering 
aerosol would have brightened the perturbed atmosphere by a few percent, and In- 
crease which would be easy to detect. Similarly, a unit optical depth of an 
aerosol with single scattering albedo .95 would have caused a decrease larger than 
any changes observable in the Figure 1 data. If an aerosol with a single 
scattering albedo very close to .984 was present, the reflectivity of the plume 
would be very nearly Independent of optical depth. So a wide range of optical 
depth of this model aerosol would match the reflectivity observations of the 
plume at 339.8 and, due to their s1milarity,331 .2nm. That is, the difference 
caused by the plume was very small. 

Figure 2 also depicts the effect of aerosol on the reflectivity of the 
perturbed atmosphere at 380. Onm The maximum aerosol optical depth consistant 
with the 339.8 and 380. Onm data is about 1.50 where the aerosol has a single 
scattering albedo of about ,984. Aerosol of this optical depth with a different 
single scattering albedo would change the reflectivity at 339.8 nm in a manner 
not observed. On the otherhand, given this single scattering albedo, the aerosol 
optical depth is limited by the observed increase in the albedo at 380. Onm, 

The effect of aerosol upon the reflectivity at the shorter wavelengths which 
are strongly absorbed by S02 is now examined. The reflectivity of the model 
perturbed atmosphere is computed using Krueger's estimates of 67matm-cm to 
compute the reflectivity at 317.5nm and 61matm-cm to compute the model reflectivity 
at 312. 5nm with no aerosol present. Different amounts of SO? were needed to 
explain the observed reflectivity for the 339.8-317. 5nm and 331 .2-312. 5nm pairs. 

This discrepancy has been attributed primarily to uncertainties in SO 2 absorption 
coefficients. When the model aerosol described above is added keeping Krueger's 
estimated amounts of SO 2 for each of the short wavelength reflectivity computation, 
the reflectivity is observed to decrease in both cases as shown in Figure 3. By 
decreasing the amount of SO 2 from 67 to 59.5 matm-cm and from 61 to 59.5matm-cm 
the reflectivities computed at 317.5 and 312. 5nm, respectively, can be matched 
with a volcanically perturbed stratum which is consistant in its properties for 
all wavelengths. This layer would contain 59.5matm-cm of SO 2 , and a model 
aerosol with an optical depth of 1.5, a single scattering albedo of .984, and 
an asymmetry factor of .83. 

The primary result of this work is a demonstration that the errors introduced 
into Krueger's TOMS inference of SO 2 amountare in the nieghborhood of 10% or less 
and not a factor of 3 or 4 which one might expect in light of other investigations 
which dealt with measurements of, or indirect inferences of, the amount of strato- 
spheric SO 2 observed after the El Chichon eruption. (see Capone, et al.(1983) or 
Evans and Kerr(1983)) 

Another interesting result is that i neons i stand es in measured SOp amount 
could result from aerosol effects;’ however, this does not rule out inaccuracies 
in S0 2 coefficicients as the source of inconsistency in this case. The aerosol 
optical depth and single scattering albedo derived in this investigation are 
consistant with the observations of optical depth at Hawaii reported by 
Del ui si et al . ( 1 983) . They report an optical depth of about .75 at 425nm 
at that date, a time when the cloud was probably somewhat dispersed. This 
compares favorably with the maximum depth of 1.5 derived in this study for 
the April 5 case. 
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SCAN NUMSEH 

Figure 1. Effects of El Chichon cloud on TOMS reflectivities. 
Logarithm of spectral reflectivity along the same trace for 
five wavelengths. (Adapted from Krueger f 1983) The dotted line 
interpolating the unperturbed reflectivity at 380nn was added 
by the author. 




Figure 2. The effect of plume aerosol optical depth and sinqle 
scattering albedo on computed atmospheric reflectivity at two 
TOMS wavelengths. The arrow Indicates the aerosol optical depth 
at which a further increase in the aerosol loading would cause 
the computed aerosol loading to exceed the observed reflectivity 
Increase at 380nm. J 


COLUMNAR S0 2 (motm«cm) 

Figure 3. The dependence of atmospheric reflectivity on columnar 
SOo The plume layer is assumed to contain an optical depth of 
1.5 of aerosol with a single scattering albedo of .984. Horizontal 
dashed lines show computed model reflectivities without aerosol 
at 312.5 and 317. 5nm with 61 adn 67 matm-cm of $0 2 , respectively. 
The plotted points show the effect of variations in columnar S02, 
when the model aerosol is present. 
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Nonlinear Analysis of the H1MAT Outer Wing 

by Larry Nelson 

Cal Poly State University 
Professor, Mechanical Engineering 


The H1MAT outer wing has, as the primary load carrying structure, 
aeroelasti^lly tailored, laminated graphite/epoxy skins. Lamlnatt coupon 
test data s»jw that the wing has significant untapped load carrying 
potential. To utilize this potential requires. In part, that a nonlinear 
structural analysis be performed. 

A preliminary, nonlinear analysis has been performed this summer. The 
basis for this analysis is an analytical representation of the experimental 
stress-strain behavior of the ply material— In this case Hercules 
AS/3501-5. The analytical representation Is In the form of a Ramberg- 
Osgood lit (references 1 and 2) to the experimental data. The analytical 
curves are differentiated to give local tangent moduli, allowing the 
nonlinear problem to be replaced by a series of piece-wise linear problems, 
each of which may be solved by classical laminate plate theory, then 
appropriately added together to give an approximation of the nonlinear 
behavior of the laminate under load (Re erence 3). Figure 1 compares the 
results of both a linear and nonlinear analysis with coupon test data 
acquired on the HiMAT laminate, (50/-50/35/-50/50) 2S . It is seen that the 
nonlinear analysis predicts quite well the coupon behavior. Based on the 
same failure criteria of maximum shear stress equal to 10 ksi, the 
nonlinear theory also better predicts laminate failure than does the linear 
theory. 

The above desor'ioed analysis can only be as good as the data base; 
i.e., the experimental stress-strain behavior of the individual plys. In 
the case of Figure 1, only the nonlinear ply shear behavior was estimated 
from (+45) laminate data of Reference 4. It is believed that significant 
nonlinear ply behavior also exists in transverse compression. At the 
present time definitive, reasonably unambiguous stress-strain data for the 
ply material have not been located. A literature search continues. 


The results of Figure 1 give confidence that a nonlinear analysis of a 
complex structure such as the HiMAT outer wing can be carried out. The 
baste approach will be to perform a series of piece-wise linear NASTRAN, 
Incremental loads analyses. In which the stress-state of perhaps 100-150 
locations In the composite wing skins will be monitored as progressively 
larger loads are (analytically) applied. The properties of the composite 


material at each of these locations will be modified in accordance with the 
current stresses. It Is anticipated that this effort will be continued by 
the writer upon his return to Cal Poly In a joint effort with Ames Dryden, 
Aerostructures Branch. It Is foreseen that the H1MAT outer wing will be 
tested In the Ames Dryden Loads Test Facility, allowing a comparison 
between test and theory. 
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IN VIVO CHANGES IN NMR PARAMETERS OF RAT MUSCLE WITH EXERCISE 


Felix S , Palubinskas, Ph.D., M.D. 

10 August 1984 

The object of the present investigation was to find out if 
measuring the NMR parameters (Tl, T2, and H) of the fore limb 
and hind limb muscles of rats, in vivo, can tell us something 
about how these muscles respond to different exercise regimens. 

In particular, we wanted to find out if the suspended rat model, 
which appeared to simulate the condition of weightlessness for 
bone, would also simulate the condition of weightlessness for 
muscle. 

A number of model systems have been used in an attempt to 
simulate weightlessness (1). The reasons for using land-based 
model systems are the following: they reduce the expense and 

difficulty of performing experiments, and they increase the 
frequency and convenience of experimentation. 

The Soviets have found that muscle atrophy occurs in rats 
in space flight and in ground-based, immobilized rats (2,3). 
Americans have found that the changes found in bone in a 
suspended rat model appear to simulate those found in space 
flight (4). 

M ATERIALS AND METHODS 

In our experiment we studied three sets of rats by proton NMR: 

(1) rats whose hind limbs were suspended , 

( 2 ) control rats, and 

(3) rats run in exercise cages (we called these rats runners ) . 

By including the set of runners, we widened our investigation 
beyond that of just the determination of the effect of hind limb 
suspension to include the effect of exercise history. 

The explanation of the details of how an NMR image is obtained 
using only radio waves incident on an animal placed in a strong 
magnetic field is given in Davis et al (5), and will not be 
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discussed here. Suffice it to say that this is an established 
technique for obtaining images and NMR parameter values of a 
human or animal body without using x-rays, or any invasive 
procedure, and without harming the animal in any way. The 
benefits of this technique are substantial, as the animal subject 
can be followed serially, in subsequent or sequential studies, 
with the prospect of using the subject as its own control. 

The experiment was run for 13 days at Ames. Three sets of 
rats were individually caged: four controls, four suspended 

rats, and five runners. Two rats from each of the three sets 
were selected for imaging. The initial weights of these rats 
ranged from 138 to 156 grams. After the second day of the 
experiment, the rats gained weight almost linearly, ending uo 
with weights that ranged from 210 to 257 grams on the 13th day. 
Suspended rats ended up heavier than the controls which in turn 
ended up heavier than the runners. 

After completing thirteen days in the experiment, the rats 
were transported from Ames to the Radiation Imaging Laboratory 
of UCSF at Oyster Point and imaged. 

It took about an hour to image each rat. For each body 
region imaged, a five-slice sequence was taken using two 
repetition times (0.5 and 1.0 sec) for each of which two spin 
echos were recorded (at 28 and 56 millisec). Four averages were 
made for each of the four possible conditions. Subsequently, 
each of the images was called up on the computer and a region 
of interest enclosed using a trakball. The computer then 
automatically computed the NMR parameters. 

R ESULTS 

The NMR data were plotted in the T1-T2, Tl-H, and T2-H planes. 
When the data were so mapped, any differences that existed between 
the corresponding muscle regions of the three sets of rats 
showed up clearly. In addition, the NMR values for ibe fore 
limb and hind limb muscles of each individual rat could be 
compared, permitting the rat to serve as its own control. 
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Examination of the maps showed: 

In the T1-T2 plane : 

(a) The hind limb muscles of suspended rats clearly map in a 
different region of the plane than the hind limb muscles of 
controls and runners . 

(b) The fore limb muscles of suspend ed and control rats map 
roughly in the same part of the plane, whereas the fore limb 
muscles of runners map in a distinctly different part of the plane. 

In the Tl-H plane ; 

(a) The hind limb muscles of suspend ed , control , and runner rats 
map in widely separated regions of the plane. 

(b) The fore limb muscles of suspended and control rats map in 
roughly the same region of the plane, whereas those of the runners 
map in a widely separated region of the plane. 

(c) The fore limb muscles of controls , and particularly the 
ru nners , map in substantially different: regions of the plane 
than their corresponding hind limb muscles. The differences 

for the fore limb and hind limb muscles of suspended rats are slight. 

In the T2-H plane : 

(a) The fore limb muscles of all three sets of rats map in 
different regions of the plane than their corresponding hind 
limb muscles. 

(b) The forelimb muscles of r unners map in a substantially 
different part of the plane from those of controls and suspended 
rats. 

DISCUSSION ' 

What does all this mean? 

It means that the exercise state of rats is reflected by 
their NMR parameters. Yet, although the measurements of NMR 
parameters are quantitative measures of the behavior of the 
so-called mobile protons in the plane of the tissue imaged, 
it is not at all clear today what precisely this means in terms 
of the molecular structure of the tissue. Nevertheless, it is 
our impression that by following the change of the NMR parameters 
of the muscles of an individual animal, the exercise state of 
various muscle groups, their development, or their state of 
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"exercise fitness" can be measured quantitatively, in vivo, 
and with complete safety. 
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To establish that the suspended rat model provides, for 
muscles, a satisfactory analog of the conditions of weightlessness 
in space flight will require NMR experiments involving space 
flight. Toward that end, an Intention To Propose document has 
been submitted to NASA Headquarters for such experimentation. 

In addition, a paper describing the present investigation is 
now being prepared and will be submitted for publication. 
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He Refrigeration for Space Applications 


Andres F. Rodriguez 


Professor, Physics Department 
University of the Pacific 
Stockton, CA 95211 


s 


The low temperature provided by 3 3 He refrigerators (0.3 - 3°K) have useful 
space applications. However, while J He refrigeration has been used on Earth for a 
number of years and in aircraft and spin-stabilized rockets, it has yet to be used 
in the low gravity environment of space. 
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The low temperature and the low surface tension of He require special design 
considerations for a space-compatible J He refrigerator. These considerations 
include: 

3 

a) The rfeed for small pores in a metal matrix to condense He vapor and to 
contain the liquid (capillary confinement). 

b) The effects of bubble nucleation and growth. 

c) The effects of the thermal conductivity within the matrix. 

These design considerations, along with a possible confinement system have been 
studies since last summer by this fellow and his NASA-Ames colleague Dr. Peter 
Kittel. 

During last summer, work was done on the validity of the model of the porous 
metal matrix proposed by Robertson et al . Also, the investigations of Ostermeier 
and his group on capillary confinement of cryogens in a metal matrix were used as 
the basis for this research project. The details of this work were included in 
the 1983 Summer Report. 

As an outgrowth of the first year summer NASA-Stanford Program, this fellow 
received a NASA-Ames Consortium Agreement Grant. This grant provided support for 
the fellow to continue his work at Ames (September - December 1983) with Peter 
Kittel while being on sabatical from the home institution. 
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The work on design considerations for a He refrigerator for space applica- 
tions was continued at the house institution during the spring of 1984 with 
frequent trips to Ames for consultation and finished during this summer. 

3 

It has been concluded that it is feasible to design a He refrigerator that 
will operate in the low gravity environment of space where accelerations can occur 
in random directions. The key to the design is to fill the liquid chamber with a 
high-conductivity porous material such as sintered copper. The pores allow 
surface tension forces to contain the liquid while the conductivity ensures that 
the evaporation occurs at the surface of the matrix rather than internally. The 
matrix also provides a favorable place for condensation to occur and suppresses 
bubble nucleation and movement. 


This research project includes mathematical and physical analysis of the role 
and strength of surface tension forces, use of the Marangoni Effect, effect of 
lateral accelerations described in terms of the Bond number and prevention of 
bubble nucleation and motion of bubbles. Finally, a discussion on how these 
considerations can be applied to the design of a refrigerator are included. One 
paper was presented at the National Conference of the ASEE and published in the 
Proceedings (3, 1332-1334, 1984). In addition, a technical memo was published by 
NASA (TM-8597T, 1984) and a paper is being prepared for publication. 

Dr. Kittel is attending an international conference in Berlin where he will 
present some of our results. 

Additional work is being done this summer on the pumping speed of cryogens 
through the pores of a matrix. Also, we are addressing other problems related to 
the surface tension of J He - ^He mixtures and the presence of surface levels in 
solutions of J He in 4 He. 

A detailed treatment of our work could be found in: 

1) Rodriguez, Andres F., "Capillary Confinement of Cryogens in a Porous 
Metal Matrix", Proceedings 1984 Annual Conference ASEE , 3_, 1332-1334, (1984). 

- 2) Kittel, Peter and Rodriguez, Andres F., "Design Considerations for a 

J He Refrigerator for Space Applications", NASA TM-85973 (1984). 

From a professional point of view, this program has been very rewarding for 
this fellow both for the research work that will be continued and for the impact 
that it will have in my teaching performance. 
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Analytic Properties of Upwind Schemes 
for The Euler Equations 

Charles W. Schelin 


The Euler equations amount to a system of conservation laws that 
model the inviscid flow of a perfect fluid. Numerous iterative schemes 
have been devised to solve discretized versions of these equations. 

Upwind schemes offer stable solutions to the difference equations by 
separating various signal directions. These schemes are usually classi- 
fied as either flux-splitting or flux-difference splitting, depending 
on the approach used to separate the signal directions. 

A number of desirable properties for a flux-splitting scheme have 
been given by van Leer [3]. Among these properties is a degeneracy condi- 
tion designed to resolve stationary shock structures in a manner similar 
to the flux-difference scheme of Osher [2]. Compering these schemes I 
found both resolve the numerical Riemann problem: 

u j = u L for j < j 0 , 

u j = u n for j > j 0 , 

L «• R — — 

where u > u > u , u being a sonic point, by yielding u. = j. 

Jo 

Using ideas similar to those developed during my first summer in 
the ASEE program, I devised a nonlinear relaxation scheme for a first 
order discretization of the one-dimensional Euler equations using the 
van Leer splitting. At some future time I hope to generalize this scheme 
to second order discretizations in two dimensions and general coordinates. 

Dennis Jespersen of NASA-Ames has shown that three of van Leer's 
properties; interpolation, degeneracy and differentiability, are incon- 
sistent in general coordinates. A flux-splitting possessing these three 
properties exists if and only if the coordinates are orthogonal. Hoping 
to draw some conclusion about flux-difference schemes I examined the 
Osher scheme. 

Flux-difference schemes depend on a numerical flux f(u,v) os well 
as the physical flux f(w). Harten, Lax and van Leer [1] use the form 


f(u,v) = Cf(u) + f(v) - d(u,v) ]/2 , 
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where d(u,v) represents an absolute flux difference, to define the num- 
erical flux. 6/ splitting d(u,v) into 


d(u,v) = d + (u,v) + d”(u,v), 


where 

d + (u,v) = f ( v) - f (u, v) 
and 

d~(u,v) s f (u) - f (u,v) , 


I was able to define the local splitting 


f + (v) = f(u*) + d + (u*,v) 
f"(v) = - d”(u*, v) 


for a sonic reference point u*\ Using this local splitting, I was able 
to show that the Osher scheme is differentiable at u* if and only if the 
coordinates are orthogonal. 

Other flux-difference schemes can be shown to be differentiable. 
Specifically, any scheme satisfying the Harten, Lax, van Leer upwind 
condition is differentiable. 

Since certain schemes which are smooth in Cartesian coordinates 
lose this smoothness in general coordinates it seems necessary to re- 
examine the transformed coordinate approach to flow problems in general 
geometries. Perhaps metric variables and flow variables can be separated 
in a fashion that will allow schemes to retain desirable properties. I 
plan to study this question in the future. 
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Preliminary Abstract Report 

Vestibular Research Facility (VRF) - Design, Control 
and Safety Analysis of the Vestibular Research 
Centrifuge and Linear Accelerator 


William Stamets, Adjunct Professor 
University of Wisconsin - Parkside, Kenosha, Wisconsin 


As a second-year Faculty Fellow at Stanford/NASA-Ames, I performed the 
following space-related research alignments during the summer 1984: 

1. Preliminary Safety and Hazard Analysis Report on Vestibular 
Research Facility (VRF). 

2. Stress Analysis Review of Impact, Seismic and Dynamic Loading 
of Experimental Centrifuge. 

3. Collected and reviewed 21 library references pertaining to hazards, 
reliability, safety, codes and standards for industrial facilities. 

4. Provided design assistance for the Linear Sled Research Accelerator. 

Copies of the above-mentioned reports will be available in about 1 week. 

The Prel iminary Safety and Hazard Analysis Report deals with the identifi- 
cation of hazards and postulated accidents at the new VRF, Building N-242 
and sets forth means for minimizing effects. It analyzes standard industrial 
equipment (crane, pressure vessel - air receiver tank, electric controls, 
etc.) and special purpose machinery (centrifuge and linear accelerator). 

The report is required by the NASA/Ames Safety Department before the facility 
can be activated. 

The Stress Analysis studies examine potential impact of loose parts flyirg 
from the centrifuge and injuring personnel and adjacent equipment. A seismic 
study is included, which considers earthquake damage in accordance with 
"Seismic and Wind Design Criteria" published by Lawrence Livermore 
Laboratories. 

An';Otated copies of pertinent sections of the 21 Codes and Standards 
references are to be pruvided with the Hazard Analysis Report. Attention 
is drawn to those portions of the references which particularly affect the 
VRF. 

My student, Daniel M. Eggert of Kenosha, Wisconsin, is here under an NASA/ Ames 
student grant which we arranged last December. He is working on the drive 
design for the Linear Accelerator Sled andds Principal Investigator. I am 
overseeing his work. We participated in a Poster Show Friday afternoon 
July 20th at the Ames Cafeteria where performance characteristics of an air- 
bearing for the sled was demonstrated. Of some twenty exhibits, eight were 
chosen to be shown at the Dryden Flight Research Facility and Dan flew down 
there August 9th. 
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During the summer, I attended and participated in the following educational 
activities given by Stanford: 

Computer Science CS-104 - Programming in PASCAL 
Mechanical Engineering MF -219 - Robotics and Computer Vision 
Aero/Astro AA-2995 - Seminar in Review of Current Research 

(Thursday evenings) in Aeronautics, Astronautics & Life 

Sciences. 

Along with these, I participated in the Monday/Wednesday noon-time Research 
Review sessions and gave a 20-minute presentation on my work program July 
11th. Also, I attended all of the Friday noon ASEE Summer Faculty Seminars 
which described a wide variety of space-related engineering subjects. All of 
these will be helpful to me in giving a talk next Fall when I have been 
invited to speak before our local engineering society. 


Will iam Stamets, P.E. 
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Analysis of Aerosol Particles by Expert System 


Scott A. Starks 


University of Texas at Arlington 
Associate Professor of Computer Science and Engineering 


Investigations into the concentration levels of aerosol particles in 
the earth's upper atmosphere have been conducted at Ames Research Center for 
a number of years. In particular, much research has been focused upon the 
study of sulphuric acid particles. Collection of aerosol particle samples 
in the stratosphere is accomplished by use of the Ames Wire Impactor which 
is mounted beneath the wing of a NASA U2 research aircraft. Current methods 
being used for classifying, counting, sizing, and developing size distribu- 
tions of the collected samples are very labor-intensive. As a consequence, 
investigation was begun last year by the author, under the NASA/ASEE Summer 
Faculty Fellowship Program, to determine the feasibility of automating the 
current method of aerosol particle analysis. 

It was discovered that commercially available image processing systems 
did not provide enough capability to discriminate the collected aerosol 
particles from contaminants and blemishes on the sampling surface. As a 
result, an automated approach based upon artificial intelligence techniques 
was formulated by the author. Work on this automated begun during the 
summer of 1983. It was continued by the author during the Academic year 
1983-1984 under the NASA/ Ames University Consortium Agreement NCA-0R780-401. 

A system diagram of the hardware used in the system is shown in Fig. 1. 
A description of the operation of the system is summarized in the following. 
First, a sample is mounted in the ISI SX-30 scanning electron microscope 
(SEM). The SEM operator then focuses the microscope on the area of the 
surface to be examined. The analog video output signal from the SEM is then 
passed through a TV-30 signal interface. The interface converts the signal 
to a standard NTSC video format. This standard signal is then digitized by 
meansof a Quantex video digitizer. A single video frame is "grabbed" by 
the digitizer under the control of an HP 1000/2100 minicomputer. The frame 
is sampled and quantized by thfe digitizer to form a (512X512) digital re- 
presentation of the image. The minicomputer archives the data onto mag- 
netic tape. Approximately 150 images are stored on a single magnetic tape 
and are used in the analysis of a single experiment. Once all the images 
have been recorded the magnetic tape is downloaded onto a VAX 11/780 
computer. 

Resident on the VAX computer is an expert system developed under this 
fellowship program for the automatic analysis of the aerosol particles. 

This system is called APES. A data flow diagram of the software organ- 
ization of the expert system is shown in Fig. 2. The front end of the 
expert system, consisting of the filtering module, the edge detection 
module, and the image segmentation module are written in the computer 
language, C. The expert system control module is written in the computer 
language, LISP. An expert system building tool called MRS was used in the 


development of the rule based classifier. The final output from the system 
is a complete description of particle size information. A more detailed 
description of the APES hardware and software systems are presented in 

Cl], C2]. 

The major advantages of the automated approach over that of manual 
inspection are many fold. First of all, analysis by expert system yields a 
large increase in the speed of reducing aerosol particle data. Secondly, 
automated system achieve appreciably higher levels of accuracy in counting 
and sizing of particles. Lastly, the extraction of shape information by 
computer yields a more accurate assessment of particle volume. 

It has been the pleasure of this fellow to work closely with a number 
of distinguished individuals in the development of APES. This fellow wishes 
to express his debt of gratitude to Dr. Henry Lum, Ms. Libby Net! and, and 
Ms. Mary Duffy for their support and assistance in the completion of this 
project. 
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Figure 1. APES Hardware Configuration 
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APES Data Flow Diagram 


















